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SUMMARY 


A method for computing three-dimensional turbulent subsonic flow in 
curved ducts has been developed. A set of tube like coordinates is given 
for a general class of geometries applicable to subsonic diffusers. The 
geometric formulation is complex and no previous treatment of this class 
of viscous flow problems is known to the authors. The duct centerline is 
a space curve specified by piecewise polynominals with continuous third 
derivatives. A Frenet frame is located on the centerline at each location 
that a cross section is required. The cross sections are described by 
superellipses imbedded in the Frenet frame. The duct boundaries are coor- 
dinate surfaces, greatly simplifying the boundary conditions. The resulting 
coordinates are nonorthogonal . 

An approximate set of governing equations is given for viscous flows 
having a primary flow direction. The derivation is coordinate invariant and 
the resulting equations are expressed in tensor form. These equations are 
solved by an efficient alternating direction implicit (ADI ) method. This 
numerical method is found to be stable, permitting solution in difficult 
geometries using the general tensor formulation. 

Flow in the entrance region of a circiHar straight pipe was computed 
at a Reynolds number of 500. At an inlet Mach number of 0.01 the results 
compared very favorably with incompressible experimental data. For a pipe 
of non-circular corss section at a Reynolds number of 500, the inlet Mach 
number was progressively increased from 0.01 to 0.1, 0.2, and 0.3 and the 
expected effects of compressibility were observed. For the case with an 
entrance Mach number of 0.3 the calculation proceded until the flow choked. 
This behavior corresponds to that expected from solution of the initial 
value problem posed. 

A series of calculations was performed for flow through a duct which 
undergoes an S-shaped bend. These cases were r\m at a Reynolds nuraber of 
500 and an inlet Mach number of 0.1. Four different duct cross sections 
were used: circular, elliptic with shape factor of 1.5> elliptic with shape 

factor of 0.66j, and super elliptic of exponent 10.0. A highly approximate 
quasi-two-dimensional inviscid pressure formulation was used for the S-bend 
calculations. The resuDts showed that in the first part of the bend a secon- 
dary flow pattern typical of that which is to be expected was produced. A 
problem was encountered near the end of the first bend where the secondary 
flows were seen to reverse. This problem is believed to result from the use 
of the quasi-two-dimensional inviscid pressure model, thus pointing out the 
need for an accurate three-dimensional inviscid pressure field. T’^o cases 
of turbulent flow through an S-shaped bend were computed at a Reynolds number 



of 105 and an inlet Mach number of 0.1. One case had a circular cross 
section, the other an elliptic cross section with shape factor 1.1. Flow 
in an axisymmetric diffuser and in a transition duct are additional cases 
discussed for completeness. 
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INTRODUCTION 


A continuing problem in the development of the Intakes for airbreathing 
propulsion systems is the design of efficient subsonic diffusers. Not only 
is the engineer faced with building an efficient diffuser, but frequently he 
must tailor the diffuser geometry to conform to certain physical constraints 
imposed by the propulsion engine and airframe. Lacking accurate generalized 
analytical design methods, the engineer must rely almost exclusively on empir- 
ical design methods based on correlations of experimental data. In the case 
of three-dimensional inlet diffusers, the cross-sectional shape of the ducting 
must vary in the axial direction, and it is frequently necessary to Introduce 
offset bends (curved duct centerlines). The diffuser geometry is, therefore, 
complicated, but perhaps more important, the offset bends induce strong secon- 
dary flows which have important effects on diffuser performance. Because of 
the vast number of geometric and flow parameters, comprehensive experimental 
programs necessary to develop generalized correlations become very costly. 
Clearly, the availability of better analytical design tools can significantly 
reduce the time and cost required to arrive at an efficient diffuser design. 

A generalized subsonic diffuser analysis capable of being used as a design 
tool must acco\int for several physical phenomena which frequently occur in 
practical diffusers. First, the analysis must be capable of treating the 
case when the wall boundary layers are turbulent and possibly of a thickness 
comparable to the dimensions of the diffuser flow passage. Secondly, the 
analysis must account for pressure gradients transverse to the direction of 
flow which can arise because of curvature of the duct centerline. Finally, 
the analysis should be capable of treating strong interaction problems; i.e., 
problems in which the viscous flow interacts with the inviscid flow. 

Because of their complexity, and particularly the interaction which 
occurs between primary and secondary flows and viscous and inviscid regions, 
three-dimensional flows in curved ducts have been extremely difficixLt to 
analyze. Rotational inviscid flow theory has provided insight into the 
behavior of some secondary flows (Hawthorne, Ref. l), and has now been devel- 
oped to the point where solutions to the full incompressible, rotational 
inviscid equations of motion can be computed (Stuart & Hetherington, Ref. 2). 
Techniques for computing three-dimensional laminar and turbulent boundary 
layers have also been developed. Some of these are surveyed by Nash & Patel 
(Ref. 3). However, there are considerable difficulties associated with the 
synthesis of secondary flow analysis and boundary layer theory into a cohesive 
method of duct flow analysis. Not the least of these difficulties is the lack 
of applicability of three-dimensional boxmdary layer theory in corner regions, 
and the treatment of interaction between viscous and inviscid flow regions. 
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In efforts to circximvent these difficulties, Patankar & Spalding (Ref. 

4 ), Caretto, Curr , & Spalding (Ref. 5), and Briley (Ref. 6) devised numerical 
methods for solving approximate governing equations which are a more or less 
natural generalization of three-dimensional boundary layer theory. In these 
studies, solutions were computed for laminar incompressible flow in straight 
ducts with rectangular cross sections. The governing equations were solved 
by integrating in a primary flow coordinate direction while retaining viscous 
stresses in both transverse coordinate directions as opposed to only one 
direction for three-dimensional boxmdary layer theory. 

In addition to neglecting streamwise diffusion, the pressure is divided 
into a mean pressure field representing the inviscid pressure and a pressure 
correction due to viscous effects. This treatment of the pressure gradients 
permits solution of the governing equations by forward marching integration. 
Subsequently, this general approach has been used to compute incompressible 
flow in helical tubes by Patankar, Pratap & Spalding (Refs. 7 & 8). Recently, 
Ghia, Ghla, and Studerus (Ref. 9 ) have employed the numerical method of Briley 
(Ref. 6) to compute flow in a "polar" duct whose cross-sections are annular 
sectors . 

Briley & McDonald (Ref. lO) have also formulated a method for computing 
three-dimensional tiorbulent subsonic flow in curved passages, which is also 
based on a forward marching integration after neglecting streamwise diffusion. 
Governing equations were derived for flow passages whose bounding walls lie 
in coordinate planes of an orthogonal coordinate system. The numerical method 
is an adaptation of the implicit techniques developed by Briley & McDonald 
(Ref. 11 ) and McDonald & Briley (Ref. 12) for application to systems of com- 
plex nonlinear parabolic and/or hyperbolic equations. In Ref. 10, some pre- 
liminary results were given for flow in a curved rectangular duct which is 
shaped like the flow passage between adjacent blades of a turbine. The pre- 
liminary calculations were made using a rather simple turbulence model re- 
quiring an ^ priori specified mixing length. Reference 13 represents a 
further development of the method of Briley & McDonald (Ref. 10 ) and an 
extension to include a more sophisticated turbulence model based on solution 
of conservation equations for turbulence kinetic energy and turbulence dissi- 
pation. The present study represents a further generalization of this method 
to encompass general coordinates and complex tube-like geometries. Although 
a derivation of the governing equations, description of the general method 
and details of the coordinate system can be found in Ref. l4, much of this 
material is repeated here in updated form for completeness. 
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A series of cases was run, including the entrance region to a straight 
pipe at a Reynolds number of 500 and at inlet Mach numbers from 0.01 to 0.3. 

A variety of cross sections was used, including ellipses of shape factors 
from 1.0 to 2.0 and super ellipses from exponent 2.0 to 10.0. The flow in a 
duct with an S-shaped bend was calculated at a Reynolds number of 500 and an 
inlet Mach number of 0.1. A highly approximate inviscid pressure formulation 
was used for these cases. Four different duct cross sections were used: cir- 

cular, elliptic with shape factor of 1.5> elliptic with shape factor of 0 . 667 , 
and superelliptic exponent of 10.0. Two turbulent S bend cases were run with 
circular and elliptic (shape factor l.l) cross sections at Reynolds number of 
105 and inlet Mach number of 0.1. An additional series of diagnostic calcu- 
lations is also included. 
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GOVERNING EQUATIONS 


Central to the present analysis is the formulation of approximate gov- 
erning equations which can he solved hy forward marching integration in the 
direction of a "primary flow". The entire flow field can thereby be obtained 
by a sequence of essentially two-dimensional calculations. This feature of 
the method results in a substantial saving of computer time and storage com- 
pared to that which would be required for solution of the full Navier-Stokes 
equations. The equations are derived in a coordinate independent manner. A 
vector field that reasonably approximates the primary flow direction is chosen 
and then used as the basis for an approximation of the stress tensor. The 
time-averaged equations are written in general conservation law form, and then 
the approximate stress tensor is inserted to obtain the approximate equations. 
Note that this process depends only on the choice of a primary vector field, 
and not on the particular coordinate system used for the numerical solution. 
The primary vector field used here consists of the tangent vectors to a cer- 
tain family of coordinate curves that are roughly aligned with the flow geome- 
try. 


The governing equations are derived from the Navier-Stokes equations for 
compressible flow of a viscous, perfect gas. In conservation law form (Ref. 
15) and, in general curvilinear coordinates (y^, y^, y3), these equations are 
given by 


A (pAq) + -A (/?v‘yg")=o (la) 

<3t 3y ' 


for continuity and 


_d 

at 





a 

ay j 


(pv‘ V j + V ']) 


ax^ 
ay ' 




= 0 


(lb) 


for momentum. Constant total temperature is assumed, and thus an energy 
equation is not required. We have used (x^, x^, x3) for fixed cartesian 
coordinates, p for density, v = v*^ ej^- for velocity, g = det (sij) = 

I det (x^V) for the metrical determinant, and xij for the components of 

O-yJ I ^ ^ ’ 

the stress tensor in the basis ej^ @ ej . In terms of the metric, the com- 
ponents of the stress tensor are given by 
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where 



H- 

Re 





(2h) 


and 



(2c) 


for viscosity y , inverse metric Kronecker deltas 6^ 

Christoffel symbols 



and 


k I <^gim ^ ) 

2 layi dy‘ 


From the ideal gas law and the constant total temperature assumption, the 
perfect gas relation has the form 


p = ABj • P + ABq- P ■ Qij ■ V ' V j 


( 3 ) 


where ABiji and ABq are constants. In all of the above, the Einstein siimmation 
convention is assumed. That is, matching upper and lower indices are to be 
summed from 1 to 3- 

To account for turbulent transport processes, the governing equations 
are time-averaged in the usual manner for tirrbulent flows (e.g., Hinze , Ref. 
l6). This process of averaging produces txirbulent correlations which are 
conventionally termed Reynolds stresses. Certain components of viscous stress 



are removed from the time-averaged equations. The process of viscous ap- 
proximation is based upon the assumption that a primary flow direction 
exists. The present approach can be regarded as a natural extension of 
three-dimensional boundary layer theory. Unlike conventional boundary layer 
theory, however, the approximate equations are to be applicable in the in- 
viscid flow region as well as the viscous region and, thus, no approximations 
are made for inviscid terms other than those to be used for the pressure 
field in subsonic flow. A detailed account of the viscous approximation is 
given in Ref. ih. 

For subsonic flow the inviscid flow region is known to be governed by 
equations which are elliptic; that is, by equations which require downstream 
boundary conditiona for solution. In this circumstance, solution by forward 
marching integration is not appropriate, at least not without some sort of 
iterative procedirre to satisfy the downstream boundary conditions. To cir- 
cumvent this problem for subsonic flows, it is therefore assimied that the 
pressure field appropriate for irrotational inviscid flow through the passage 
represents a given, reasonable first approximation to the actual pressure 
field. Thus, inviscid axial pressure gradients computed with appropriate 
downstream boundary conditions are "imposed" upon the flow, much as in con- 
ventional boundary layer theory, so as to permit solution by forward marching 
integration for subsonic flows. For Internal flows, the inviscid pressure 
gradients are corrected for internal flow losses associated with the well 
known viscous pressure drop and blockage effects by a process which is con- 
sistent with forward marching integration. The imposition of inviscid 
pressure gradients incorporates ^ priori the elliptic effects associated 
with a subsonic pressure field without the necessity of solving elliptic 
equations other than for an inviscid flow. The inviscid pressure field can 
be generated from any convenient source and does not necessarily require the 
solution of compressible or even three-dimensional inviscid flow equations. 

In the S-shaped bend cases presented, the three-dimensional inviscid pressure 
field was not available, so a highly approximate method was used to provide 
the inviscid pressures. 


Turbulence Model 

The mixing length model used in this analysis employs the eddy-viscosity 
formulation for the Reynolds stresses, i.e.. 




V ' vJ 


d\j i 

dXj 


( 4 ) 
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Hence, this formulation still suffers from the physical shortcoming that 
there is zero Reynolds stress wherever the velocity gradient is zero. In 
addition, the eddy viscosity formulation is isotropic which may he incorrect 
in many three-dimensional and swirling flows. However, for practical cal- 
ciilations of turbulent Internal flows there are no other available transport 
models which are as suitable or even as relatively well developed. 

The mixing length turbulence model employed in this analysis is based on 
a mixing length distribution. The math^iatical form of the expression for 
the turbulent viscosity follows from Ref. IT: 

< 5 ) 


where e is the mean flow rate of strain tensor 

e= 1/2 [(Vv) +(Vvf] (6) 


The mixing length Si is determined from the empirical relationship of McDonald 
& Camarata (Ref. l8 ) for equilibrium turbulent boundary layers which can be 
written 


(y) = 0.09 Sj, [/c y / (0.09 S J • 


(Ta) 


where 6-jj is the local boundary layer thickness, k is the von Karman constant, 
taken as 0.43, y is distance from the wall, and ^ is a sublayer damping factor 
defined by 


(y'*'- y*)/ ^ 


(Tb) 


where P is the normal probability f-unction, y'*’ = y (t/p )l/2/(y/p ) , x is local 
shear stress, y"*" = 23, and = 8. 
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As a means of simplifying the present preliminary calculations, y is 
taken as the distance to the nearest wall, and 6 -^ is simply a rough estimate 
of the average boundary layer thickness in the duct. The shear stress t 
appearing in y was computed assuming a local skin friction coefficient 
of about 0.00359 the value appropriate for the inlet boundary layers. Here, 
is based on the local inviscid velodity Uj, i.e., = x/l2pUj^. 
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CURVILINEAE COORDINATE SYSTEM 


An overview of the nonorthogonal tube-like coordinate system is presented 
here. A complete description of the construction and component parts of these 
coordinates can he found in Ref. l4. 

Tuhe-like coordinates have been constructed to provide a natural setting 
for the study of flows within, between, or outside of a set of prescribed 
tubes. The prescribed boimdary tubes become coordinate surfaces, and, as a 
result, the specification of fluid dynamic boiandary conditions is greatly 
simplified. Although the equations of motion contain more terms than for a 
cartesian system, this does not add excessively to the run time of the pro- 
gram. The basic geometry of the bounding tubes then provides the intrinsic 
constraints upon the coordinate construction. Since the primary goal is the 
computation of fluid flows within nontrivial geometries and not the develop- 
ment of coordinate systems per se, the coordinates are kept as simple as 
possible, given the desired generality. 

The first step in the construction of tube-like coordinates is to create 
a suitable family of two-dimensional sxirfaces which, in some sense, are 
transverse to a given centerline. If the transverse surfaces are selected 
to be two-dimensional- planes , then the construction of coordinates is greatly 
simplified. The fluid dynamic computation is only marginally different due 
to the resulting coordinate nonorthogonality. Consequently the coordinate 
system that we shall construct will have planar transverse surfaces. 

Since each planar transverse surface is a linear subspace of the real 
three-dimensional Euclidian vector space r 3, any such plane can be completely 
specified by any two spanning linearly independent vectors in R^. The speci- 
fication of the planar family of transverse surfaces is then a result, of a 
construction of two vector fields along a given centerline curve in r 3. The 
origin of each plane is chosen to coincide with the associated centerline 
point (Fig. l). To assure that the planes are always Indeed transverse, it 
will be assumed that they are orthogonal to the centerline at their origins. 
Tube- like surfaces are generated by loops about the planar origins which deform 
in some way as we move along the centerline curve; in general, these tube-like 
surfaces will not intersect the transverse planes orthogonally. Thus, only 
the centerline direction determines the transverse nature of the cross sec- 
tional planes. Specifically, the centerline tangent vectors form a vector 
field which, at each point, is orthogonal to the plane of the two transverse 
vectors, and thus each centerline point carries a triple of linearly inde- 
pendent vectors. By the Gram-Schmidt orghogonalization procedure, each such 
triple of vectors can be made into an orthogonal set, and hence, an ortho- 
normal set which is simply called a frame. Thus, tube-like coordinate systems 
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Ip ^ 

notation, y = 0 , y =r and y -t for pseudo-angular, pseudo-radial, and axial 
variables. In this notation, we have thus far developed (l) a length factor, 
L=L(y^, y^, y^), which is a generalization of radius, (2) an angular distri- 
bution function, 0=0 (y^, y^) which is a generalization of angle,^(3) a^ 
r ot at ion_^f unction, (y3), and (l^) the Frenet frame, (Vp, V2 j ^ 2 '^= (Vj(y^), 

V2(y^)> V^(y3)) upon which the coordinates are built. That the length factor, 
L, and the angular distribution function, 0, give us a generalization of polar 
coordinates is obvious since polar coordinates are easily retrieved by taking 
L(y^> y^5 y^) = y^ and ©(yl, y3) = yl. It is also worth noting that the 
angular distribution function, 0, was chosen to be independent of pseudo- 

P 

radius, y . Although it is not immediately evident, we have removed a 
considerable amount of potential computational complexity in the process of 
obtaining metric information by limiting the number of derivatives which must 
be computed. Furthermore, there is no real loss of flexibility in the con- 
struction of angular distribution functions. Since most commonly used 
analytic descriptions of loops are, in fact, controlled by a collection of 
parameters which depend only on axial location, y3, a knowledge of only these 
parameters is often sufficient for the construction of the angular distribu- 
tion function. For example, if the loops were to consist of a family of con- 
centric homogeneous ellipses, then the major and minor axes of the outermost 
ellipse would form a collection of two such parameters. 

With the above functions and the Frenet frame, the class of tube-like 
coordinates comes directly out of the transformation 

X = 7 + L cost/) + V 3 sin </) I (8a) 


V “I Q *5 

which transforms curvilinear coordinates, y = (y^, y'^, y^) into cartesian 
coordinates x = (x^, x^, x^) where 


</) (y', y 3) = 0 (yi ,y 3) + lUy3) 


(8b) 


At each transverse location, y8, the space curve vector, Yj translates the 
origin to the space c;orve. At a given pseudo-angle, y^, a unit vector, V2 
cost)) + V3 sin<f>, is determined by the sum, cj) = 0 + of the radial distribution 
function, 0, and the transverse rotation, fi. This unit vector sweeps out a 
full 360 degs in the transverse plane as y^ passes through all of its values. 
Hence, we could call this a direction pointer for the transverse plane. When 
this direction pointer is scaled by the length factor, L, we obtain a point 
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are constructed from a specified centerline curve and an associated frame 
field. Now the basic question is whether there is a canonical construction 
of tube-like coordinate systems from either a given centerline or a given 
frame field. From the theory of space c^lrves (Ref. 19 ), it is well known 
that for positive ctirvature and specified torsion there is a local one-to-one 
correspondence between frame fields and space curves which pass through a 
given point. Thus, for nonzero curvature, the centerline space curve has a 
canonical frame field which is known as the Frenet frame. Consequently, the 
coordinates will be derived from the Frenet frame when it exists. At center- 
line points of zero ciirvature, the Frenet frame is degenerate and must be 
treated specially. 

Once the Frenet frame of the_^ space jcurve y has been established, the 
unit normal and binormal vectors V2 and at each point of y determine a 
transverse plane orthogonal to the unit tangent vector (Fig. 2 ). Relative 
to any such transverse plane, these vectors are also the standard orthonormal 
basis. Consequently, we can examine the plane separately from the curve, y, 
which will only appear as the point at the origin. In two dimensional func- 
tional terminology, the unit normal direction can be considered as the abscissa 
and the unit binormal as the ordinate; or more simply, as x and y axes, re- 
spectively. Since the tube-like coordinates are to be generated from some 
family of tubes encasing the space ciirve, y , a cross-sectional cut by a 
transverse plane produces within the plane a family of loops about the origin. 
We shall assume that each loop is representable by a strictly monotone radial 
function of angle. In this regard, a polar type of description is the most 
suitable. But, of course, the loops are usually more complicated than 
circles, and thus, we must replace the radius by a function L of both radial 
and angular variables r and 0 . Furthermore, when noncircular loops bound a 
cross section of fluid, there are regions of varying wall curvature. In a 
numerical solution, it is desirable to put proportionately more mesh points 
in regions of higher curvatiure than in regions of less curvature. Conse- 
quently, an angular distribution function, 0, is a good replacement for the 
simple angular specification 0, of simple polar coordinates. The net result 
is a generalization where polar coordinates are replaced by a pseudo-radius, 
r, and pseudo-angle, 9. Since the loops generally vary from transverse 
plane to transverse plane, the pseudo-radii and angles must also be fiinctions 
of axial location, t, on the centerline space cirrve, y{t). Since the normal 
and binormal directions are usually fimctions only of the centerline cirrve, 
y(t), our loops may have symmetries that do not reflect about either of these 
Frenet directions. Since the use of known symmetries is a great simplifica- 
tion in most problems, we need an option which allows one to define axes that 
can be aligned in an optimal way. This option is easily established from the 
specification of a function, J^(t), which is a rigid rotation relative to the 
normal-bi normal directions. To bring this development of tube-like coordinates 
within the framework of the preceding tensor derivations, we shall use the 


13 


of our transformation. Since the length factor depends on all three varia- 
bles, any set of tube-like surfaces can be obtained provided, of course, that 
loops are representable by a strictly monotone radial function of angle and 
also that no two transverse cross sections are allowed to intersect. 

In a geometric setting, the transformation is really an embedding of 
tube-like coordinate systems into three dimensional Euclidian space. An 
illustration is provided in Fig. 3. From the transformation, it is also easy 
to see that the surfaces of constant y3 are the transverse planes , the sur- 
faces of constant pseudo-angle, y^, are ruled surfaces generated from the 
centerline curve, y, and the surfaces of constant pseudo-radius, y^, are 
just the concentric tubes about the space curve, Y- Separate illustrations 
of these various coordinate surfaces are given in Figs. 4a, 4b, and 4c, 
respectively. 


l4 


I II 


nil II 


METHOD OF SOLUTION 


Background 

The governing equations can he solved (after modeling the Reynolds 
stresses) following the general approach developed hy McDonald & Briley (Ref. 
12) for laminar supersonic flow in rectangular jets. A detailed discussion 
of the calculation procedure is not included here, as such a discussion would 
be lengthy, and discussions of the general approach are available elsewhere 
(Briley & McDonald, Ref. 11; McDonald & Briley, Ref. 12). The method used is 
based on an implicit scheme which is potentially stable for large step sizes. 
Thus, as a practical matter, stability restrictions which limit the axial 
step size relative to the transverse mesh spacing and which become prohibitive 
for even locally refined meshes (e.g. , in laminar sublayers) are not a factor 
in making the calculations. The general approach is to employ an implicit 
difference formulation and to linearize the implicit equations by expansion 
about the solution at the most recent axial location. Terms in the difference 
equations are then grouped by coordinate direction and one of the available 
alternating-direction implicity (ADI) or splitting techniques is used to 
reduce the multidimensional difference eqiiations to a sequence of one-dimen- 
sional equations. These linear one-dimensional difference equations can be 
written in block-tridiagonal or a closely related matrix form and solved 
efficiently and without iteration by standard block elimination techniques. 

The general solution proced\are is quite flexible in matters of detail such 
as the type and order accuracy of the difference approximations and the 
particular scheme for splitting m\ilti dimensional difference approximations. 
Based on previous experience of the authors, however, it is believed that 
the consistent use of a formal linearization procedure, which incidentally 
requires the solution of coupled difference equations in most instances, is 
a major factor in realizing the potential favorable stability properties 
generally attributed to implicit difference schemes. 


Computation of Mean Pressure Drop 

To compute the mean press\ire drop, the streamwise pressure gradient term 
in the momentiom equations is replaced by (dpjn/dx)’^ plus a function of two 
parameters having the form 


7T, (X,Z) , 7T2(x,y) 


(9) 
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where and t 72 are introduced as a computational device to allow an im- 
plicit splitting of the change in the pressure drop term across a step. The 
dependence of iri on z and tt 2 on y is removed (to order Ax) as part of the 
solution procedure. During a step of the solution procedure, the mean 
pressure drop term is replaced by 


(dpm/dx)^/2 + 7r*(x,z) - 7T2'^(x,y) 


(lOa) 


dirring the first ADI sweep and by 


(dprp/dx)^/2 - TT,'^ (x,z) + 7T**(x,y) 


(lOb) 


during the second ADI sweep. The variables tt* and tt** are parameters computed 
as part of the implicit solution procedure so as to satisfy an integral mass 
flux condition. The parameters and tt 2 contain information concerning the 

mean pressure drop for a given step. To extract this information, it is nec- 

^ . . 

essary to define mean values tti and tt 2 9 averaged in the transverse directions 
such that 



(lla) 


^2 




(lib) 


It will be seen that, at the beginning of a step, and are adjusted to 

have zero mean for the step about to be taken, i.e., ttJ = = 0. After a 

step is complete, the updated mean pressure drop is obtained from the area 
averaged parameters as 


, n + i 

(dpm/c^x) 


n 

(dpn^/dx) + 


7T| + TTg 


( 12 ) 
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and the unaveraged parameters are adjusted for the next step so as to have 
zero mean as follows : 


■w, 


n + i 


7T, - 


^ n +1 
^ 1 * ^2 


7T2 




- 7T. 




(13) 


The difference approximation across a complete step consisting of hoth ADI 
sweeps is obtained by adding Eqs. (lOa & b) and invoking the definitions (13) 
to obtain 


(dp^/dx)"'*'' 


(dpj^/dx) + 


* 4. ** 

7T, + TTg 


r n+i 
+ [tt, - 



(lU) 


The last two bracketed terms in Eq. (l^) represent error terms of order Ax 
and have zero mean. Thus, the mean values of the parameters irq and Trp account 
for the change in mean pressure drop across a step, whereas the transverse 
variation of and vp from the mean is effectively subtracted out during each 
step to within the aforementioned error terms which in practice are found to 
be quite small. The overall technique is designed to permit the implicit 
computation of mean pressure drop using an ADI method without the necessity 
of iterating. Calculations which serve to verify the procedure are presented 
subsequently. 


Integral Mass Flux Condition 

The value of vq* and Tr2** are determined as part of the splitting process 
so as to satisfy the integral continuity condition for internal flow with no 
mass addition 


J ^ dyou/dx dydz=0 (15) 

yi Z| 


The integral condition (15) cannot be split conveniently for use with an ADI 
method; however, the presently used split form of the continuity equation is 
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(you)*- (/3 u)^]/Ax + d{pvf/dy=0 


(l6a) 


[(^u) ** - {/ou)*j/Ax + a(yow)** /dz=0 


(i6b) 


Integrating Eqs. (l6) with respect to the implicit direction and invoking the 
condition of zero mass addition yields 


[(^u)*-(^u)"]dy = 


0 


(17a) 


f ^ [(/3u) ~(/>u) j dz = 0 


(1Tb) 


which is the split form used as auxiliary conditions to determine and tt 2 - 
Integrating Eqs . (IT) with respect to the explicit direction yields 

/ [(yOu)*- (you)^l dydz =0 (l8a) 

•'z, •'y, L -I 


[(you)** “(yOu)* ] dzdy =0 (i8b) 

yi ^1 


and thus the integral continuity condition is satisfied across a complete 
step. 
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Solution of the Split Difference Equations 


The coupled set of linear implicit difference equations arising along 
rows of grid points during each step of the ADI solution procedure, together 
with the prescribed boundary condition, can be written in a form ha.ving the 
following matrix structiore 



For each grid point index i , <t>i is a column vector containing the dependent 
variables P, u, v, w. P, a single variable, is the pressure correction pq or 
P 2 depending on the ADI step. Aq , Bq , and Cq are square (i^x4) matrices 
containing the implicit difference coefficients. fq is a column vector con- 
taining the implicit coefficients of P, and dq is a column vector containing 
only computationally known quantities. There are N+1 grid points along the 
row under consideration. Difference approximations for the four governing 
equations are associated with symbols having subscripts 1 through N-1; the 
subscripts 0 and N are associated with the boundary conditions, which may 
involve up to three grid points. Equation (19) represents 4 (N+1) linear 
equations in 4 (N+l) dependent variables plus the additional parameter var- 
iable P. Thus, one additional boundary condition or auxiliary relation is 
required to close the system, and this additional relation is presumed to be 
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available when needed. Excluding the elements Cq, Ajj, and the matrix 
structure of Eq. (19) is block tridiagonal, and direct solution by standard 
block elimination techniques (cf., Isaacson & Keller, Ref. (20) is both 
straight fo 2 Tward and efficient. The precise scheme used here consisted of 
Gaussian elimination for a simple tridiagonal system (sometimes called the 
Thomas algorithm) but with elements of the tridiagonal matrix treated as 
square submatrices rather than as simple coefficients. The required in- 
verses of diagonal submatrices were obtained by a Gauss-Jordan reduction. 
The additional operations necessary to include the non-block-tridiagonal 
elements Cq, Ajj, and f^ are easily incorporated provided the original block 
tridiagonal coding is carefully organized. 
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APPLICATION TO INTERNAL SUBSONIC FLOW 


The analysis of the previous sections has "been applied to the flow in 
a pipe whose centerline is specified by piecewise polynomial functions with 
continuous third derivatives. Cross sections of the pipe in planes normal 
to this centerline are described by the equation of a superellipse (in local 
cartesian coordinates): 


b°|: 


IyI 


m 


( 20 ) 


where a is the superelliptic e^cponent constrained to the range 2< a <10, 
b is the ratio of major to minor axes (sometimes called the shape factor), 
and m is the length of the minor axis. 

Throughout the remaining discussion, all variables in the governing 
equations are nondimens ional , having been normalized by the following ref- 
erence quantities: distance, L-p ; velocity, U^; density Pi>; temperature, Tp ; 

P P 

total enthalpy, Uri pressure, PpU^; vxscosity, y^. Here, the subscript r 
denotes a reference quantity. After the various modifications outlined above 
are implemented, the governing equations (2a-e) can be written, after neg- 
lecting all density fluctuations, as 


ay' 


d 

d y j 


pv ' V ^ 


+ 


g'j p + oJ'j) 




(21b) 


where is the viscous portion of the stress tensor (Eq. 2a) modified by 
the viscous approximation of Ref. l4, a and B are redefined to include the 
Reynolds stresses as: 


a 


ij _ 

k ■ 


Re 



+ 


ayk/ 


(22a) 
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and 


Hf - - gi«S') (sat) 


and the pressure is 


P = 


r-i ^ 


P [e- 


T V v'] for j = i,2 


(23a) 


P = Pi + Pm for j = 3 


(23b) 


where y is the ratio of specific heats, Pp is the specified inviscid pressure, 
and Pjji is the mean pressure drop computed by the technique outlined above. 


Boundary and Initial Conditions 

In computing turbixlent flow near solid boundaries a problem can arise 
from the juxtaposition of large and small physical length scales within very 
small physical distances. One method of solving this problem is to use a 
coordinate transformation to pack computational points near the wall to re- 
solve the sublayer. Spacing between grid points is increased as the distance 
from the wall and the length scale increase. A second method is to employ an 
analytical wall function to model the near wall region, thus avoiding solving 
the governing equation in this region (Refs. 8 and 21 ). Both methods are 
available within the computer code. Because of the economy of grid points, 
and the resulting reduction in computational cost, only the wall function 
method has been used to date. 

To implement the wall function boxmdary conditions, the first grid point 
away from a wall is positioned in a highly turbulent region of the flow. 

Since the near wall region is not resolved, the usual specification of no- 
slip conditions is no longer satisfactory. This is shown schematically in 
Fig. 5 by reference to a typical turbulent boundary layer velocity profile 
with a superimposed grid point distribution. The line segment CD represents 
the slope of the velocity profile at the first grid point. If one-sided 
differencing were used for 3V/9n together with the no-sllp condition at the 
wall, the computed slope would be the line segment AB and substantial error 
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wolild occur. To circumvent this difficulty, a fictitious slip velocity O' 
is introduced and is determined from the wall function boundary condition 
such that the finite difference value of 3V/9n at the first grid point is 
consistent with the ass-umed analytical form for the turbulent velocity 
profile. 

To derive wall function boundary conditions, it is assumed that the 
velocity component normal to the wall is zero and that the magnitude of the 
resxiltant velocity vector ^ u, obeys a universal log law 


u=Ur j^n (urn/v)/K-+c 


(24) 


1/2 

where is the friction velocity U-^ = (xy/p ) , k = 0.4l is the von Karman 

constant, and C is a constant taken as 5-0. The boundary condition is ob- 
tained by differentiating Eq. (24) with respect to n and neglecting the 
contribution due to skewing of the resultant velocity u, yielding: 

<311 

(25) 


The wall density is specified from a knowledge of the pressure. At the 
symmetry planes of 9 = constant, symmetry conditions are used. 

Since a small tube is placed around the centerline to circumvent the 
geometric singularity there, careful handling of the boundary conditions is 
required. Since this area is distant from the walls, gradients of the fluid 
properties are expected to be small. Consequently the first derivative of 
the streamwise velocity is set to zero. The pressure gradients are all mild 
so that the second derivative of the density is also set to zero. The cross 
flows at the centerline tube must be set to insure that total mass flow is 
conserved. The transverse velocity at the centerline tube is set as a func- 
tion boundary condition whose value is equal to the average transverse velo- 
city at the previous computational plane in the neighborhood of the center- 
line. 


To obtain upstream starting conditions, the inviscid velocity distribu- 
tion in the plane at the initial axial location x = xq was scaled near the 
wall by normalized two-dimensional compressible turbulent boundary layer 
velocity profiles having a prescribed thickness, Reynolds number, and Mach 
niimber. The particular profile shape used consisted of a Coles (Ref. 22) 
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profile modified for compressibility as suggested ty Maise and McDonald 
(Ref. 23) and made continuous through the sublayer region following the 
suggestion of Waltz (Ref. 2 h) . With the backward difference formulation, 
starting conditions for v and w are not required other than in the lin- 
earization procedure, and setting v = w = 0 at the starting plane was found 
to produce a smooth and reasonable starting process for the calculation. 
Given the velocity, the density can be computed from the inviscid pressure 
Pj, here assumed constant across the Initial value surface. 



DISCUSSION OF RESULTS 


Results from computed solutions are presented for two flow problems. 

The first problem is the laminar flow in the entrance region of a straight 
pipe. This flow has been thoroughly studied both experimentally and analyt- 
ically providing a good means for testing the validity and accuracy of the 
overall calculation procedure. The second problem treated is that of flow 
through an S-shaped bend. Calculations have been performed at a series of 
Reynolds numbers, both laminar and turbulent, for several cross-sectional 
duct shapes . 

For completeness, some results from a series of diagnostic runs are 
also included in this report. These cases consist of flow in an axisym- 
metric diffuser and of flow in transition ducts whose cross sectional shape 
varies from nearly square to round. The results from these runs are not as 
yet adequate, pointing out the need for further code development addressed 
to these problems. 

The flow in the entrance region to a circular straight pipe was computed 
at a Reynolds number of 500 based on pipe radius. The initial Mach number 
was set to 0.01 to approximate the incompressible solution. The initial 
boundary layer thickness was set at 0.075 diameter. This test case is con- 
sidered important as a means for verifying large portions of the code since 
the entire nonorthogonal geometric package was used to generate coordinates 
that reduce to cylindrical coordinates for this case. The entire general 
tensor formulation of the equations was used to generate the familiar 
"parabolized" Navier-Stokes equations in the cylindrical coordinates. In 
addition the entire ADI procedure was used to split the equations into radial 
and circumferential parts for the two part solution procedure even though 
the problem was axisymmetric . The results of this calculation were compared 
to the detailed experimental data of Reshotko (Ref. 25). Excellent agreement 
between the measured and computed velocity profiles was found over the devel- 
oping region as seen in Fig. 6. In the fully developed limit, the calciLLa- 
tions, measurements and classical analytical solution were virtually identical. 

This test case was repeated for a variety of cross sectional duct geom- 
etries including ellipses up to a shape factor (ratio of major to minor axis) 
of 2.0. A series of test cases was also run for a series of ducts whose cross 
sections were superellipses with exponents ranging from 2 to 10. Although no 
experimental data were found for comparison with these solutions, the results 
nevertheless appear reasonable. Figure T presents the centerline velocity 
distributions for the series of superelliptic ducts. 
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An additional set of straight pipe calculations was run to check the 

computation of the effects of compressihility . Using a straight pipe of 

superelliptic cross section (exponent 3.0) the initial Mach number was 

progressively increased from 0.01 (incompressible solution) to 0.1, 0.2 

and 0.3. As expected, the compressible solutions did not exhibit the 

limiting value of centerline velocity. At a given nondimens ional distance 

from the start of the pipe ( ?/^ ) the nondimens ional centerline velocity 

Re 

increased as the effects of compressibility increased as seen in Fig. 8. 

This behavior is entirely consistent with the well know effect of Mach 
number on pressiu'e loss due to friction. 

A final test of the calculation procedure occurred in the case having 
an initial Mach number of 0.3. As the flow proceeded downstream, the 
pressure losses due to friction caused the centerline Mach number to 
increase past the sonic point into a region where the centerline flow 
became supersonic. The calculation remained stable in this region. Finally 
the velocity in the centerline region became sufficiently large that, in 
the mean, the flow encountered the sonic singularity. This singularity 
arises from the factor (l-M^) ^ appearing in the one-dimensional pipe flow 
equations where M is here taken to be a mean Mach number at an axial loca- 
tion. Although this term does not appear directly in the three-dimensional 
equations being solved, the phenominological correspondence between the 
computed and the one-dimensional results strongly suggests a valid integra- 
tion of the finite difference equations. When the flow in the pipe choked 
the calculation procedure broke down immediately. This behavior corresponds 
to that expected from solution of the initial value problem posed. 

The second set of results being presented is a series of calculations 
of flow through an S-shaped bend. The baseline configuration is a circular 
pipe with an offset of 1 diameter in a length of 10 diameters shown in Fig. 

9. The first test case was laminar with a Reynolds number of 500 based on 
radius, and an inlet Mach number of 0.1. The accuracy of the flow predictions 
obtained for this sequence of calculations is limited by the present lack of 
an adequate inviscid pressure field as determined by the potential flow 
through the S-bend geometry. This was not a problem in the previous straight 
duct calculations since in these geometries the inviscid pressure is known 
to be constant. 

To demonstrate the potential of the present method, a highly approximate 
method was used to generate an inviscid pressure field. Using the radius of 
curvature of the pipe centerline at each computational plane, a two-dimensional 
inviscid free vortex solution in the plane of the pipe centerline and the 
radius of curvature was assumed. The corresponding pressure at each computa- 
tional point was computed. Although this simplified model does not produce any 
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pressure gradients normal to the plane of symmetry and does not account for 
any cross flows, it does produce a pressTore higher at the outside of the tend 
than on the inside. Using the same procediare at each computational plane, 
the inviscid pressure gradients were computed from backward differences of 
these pressures. 

Using this quasi-two-dimensional pressure field the flow through the 
S-duct of Fig, 9 was computed. The results showed that in the first part of 
the bend a secondary flow pattern typical of that which is to be expected 
was produced. The low inertia fluid in the boimdary layer moves from the 
high pressxire at the outside of the bend along the wall toward the inside of 
the bend. The flow outside the boundary layer responds to the centrifugal 
forces caused by the t\irning of the duct. A series of vector plots of the 
cross flow velocities is presented in Fig. 10. A problem was encountered 
near the end of the first bend where the secondary flows were seen to reverse. 
This problem is believed to result from the use of the quasi-two-dimensional 
pressure model for this three-dimensional problem, thus pointing out the 
need for a three-dimensional inviscid pressure field. 

Several additional cases for the curved centerline in elliptic ducts of 
shape factor O .667 and 1.5 and in a superelliptic duct of exponent 10.0 were 
run to demonstrate the geometric capability of the code. In figures 10-13, 
Station No. 5 is the end of a 1 1/4 diameter long straight section. From 
Station No. 6 through Station No. 13 the radius of curvature decreases until 
it reaches 15 diameters. 

An additional series of cases was run for turbulent flow at a Reynolds 
number of 10^ and an inlet Mach number of 0.1. Using the same centerline as 
the above cases, ducts with circular and mildly elliptic (shape factor l.l) 
cross sections were r\m. The results of these calculations are not available 
for publication at this time. Each of the S-shaped bend cases was run on a 
15 X 15 computational grid at each of 45 streamwise planes. Using the 
Univac 1110 computer, each case ran approximately one hour. 

Several additional calculations were performed in the course of the 
present investigation. Although these are diagnostic calculations whose 
results are not as yet adequate, several of these are included for complete- 
ness . 


The case of flow through an axisymmetrlc diffuser produced a problem in 
which the velocity near the centerline decreased rapidly, even at wall angles 
on the order of 1°. This unrealistic deceleration was confined to the 3 or 
4 grid points near the centerline. A detailed examination of this calcula- 
tion showed that the streamwise velocity and density were changing locally 
in a compensating manner to maintain the imposed static pressure. 
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The other set of diagnostic runs is of flow in transition ducts whose 
cross section varies from nearly square to circular as shown in Fig. l 4 . 

The calculations were run both laminar and turbulent, and additional pipes 
of constant cross section were added both before and after the transition 
section. When these were being r\in, the computer program was being coupled 
to the computer graphics system developed under the direction of Bernhard An- 
derson at NASA Lewis Research Center. The lines of constant velocity at 
several cross sections produced by the computerized graphics system are shown 
in Fig. 15. Again, difficulties were apparent in the region of the center- 
line. Fiarther investigation is needed to resolve the problem encountered 
near the centerline in the axisymmetric diffuser and transition duct calcu- 
lations . 


28 



RESUME 


A method for computing three-dimensional turbulent subsonic, flow in 
ciorved ducts has been developed. A set of tube like coordinates is given 
for a general class of geometries applicable to subsonic diffusers. The 
geometric formulation is complex and no previous treatment of this class 
of viscous flow problems is known to the authors. The duct centerline is 
a space curve specified by piecewise polynominals with continuous third 
derivatives. A Frenet frame is located on the centerline at each location 
that a cross section is required. The cross sections are described by 
superellipses imbedded in the Frenet frame. The duct boundaries are coor- 
dinate surfaces, greatly simplifying the boundary conditions. The resulting 
coordinates are nonorthogonal . 

An approximate set of governing equations is given for viscous flows 
having a primary flow direction. The derivation is coordinate invariant and 
the resulting equations are expressed in tensor form. These equations are 
solved by an efficient alternating direction implicit (ADl) method. This 
mmierical method is found to be stable, permitting solution in difficult 
geometries using the general tensor formulation. 

Flow in the entrance region of a circular straight pipe was computed 
at a Reynolds number of 500. At an inlet Mach number of 0.01 the results 
compared very favorably with incompressible experimental data. For a pipe 
of non-circular corss section at a Reynolds number of 500, the inlet Mach 
number was progressively increased from 0.01 to 0.1, 0.2, and 0.3 and the 
expected effects of compressibility were observed. For the case with an 
entrance Mach number of 0.3 the calculation proceded until the flow choked. 
This behavior corresponds to that expected from solution of the initial 
value problem posed. 

A series of calculations was performed for flow through a duct which 
undergoes an S-shaped bend. These cases were run at a Reynolds number of 
500 and an inlet Mach number of 0.1. Four different duct cross sections 
were used: circular, elliptic with shape factor of 1.5) elliptic with shape 

factor of 0.66T, and superelliptic of exponent 10.0. A highly approximate 
quasi-two-dimenslonal Inviscid pressure formulation was used for the S-bend 
calculations. The results showed that in the first part of the bend a secon- 
dary flow pattern typical of that which is to be expected was produced. A 
problem was encountered near the end of the first bend where the secondary 
flows were seen to reverse. This problem is believed to result from the use 
of the quasi-two-dimensional inviscid pressvire model, thus pointing out the 
need for an accurate three-dimensional Inviscid pressure field. Two cases 
of turbulent flow through an S-shaped bend were computed at a Reynolds number 
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of lo5 and an inlet Mach number of 0.1. One case had a circular cross 
section, the other an elliptic cross section with shape factor 1.1. Flow 
in an axisymmetric diffiiser and in a transition duct are additional cases 
discussed for completeness. 
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APPENDIX 


COMPUTER CODE INPUT AND OUTPUT 


The input to the PEPSIG code is provided through three Namelist tables. 
Although the code is still in the development stage and elegant input has 
not yet been devised, the Namelist method provides a serviceable means of 
setting up the cases to be run. The three Namelists are, in order of being 
read: GEOB, LISTl, and COEFL. The nimibers in parenthesis following the 

explanation of the input variable are the recommended values. 

GEOB : 

NRAD - Number of grid points in the radial direction (l5) 

MQ - One less than the number of grid points in the circumferential 

direction (l4) 

RTUBE - Radius of the tube put around the centerline to remove the 

geometric singularity (O.O 5 ) 

VIS(2) - Fraction of a quadrant to be computed (0.5 for half quadrant 

symmetry, 2.0 for two quadrants) 

P(l,j) - Coefficients of polynominal in X for shape factor - ratio of 

major to minor axis where: 

S.F. = P(l,l) + P(l,2) X + P(l,3) + P(l,4) x3 ... 

+ P(l,ll) x^° 

P(2,j) - Coefficients of polynominal in x for radius of major axis 

P(4,j) - Coefficients of polynominal in x for superelliptic exponent 

between 2 and 10. 

PA(3,1jJ) - Coefficients of polynominal in x for traverse cartesian com- 

ponent of centerline curve. Set to zero for straight duct. 

NOTE: The P and PA arrays determine the location of the grid points. Since 

the position vector of the grid points is differentiated three times within 
the code, the piecewise polynominals describing these locations must be con- 
tinuous to three derivatives. One must be careful of these derivatives, they 
have a great effect on the code. 
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LISTl : 


YZERO - Reference length in ft (l.O) 

RZERO - Reference density in slxigs/ft^ 

UZERO - Reference velocity in ft/sec 

CMACH - Reference Mach Number 

NOTE: Set either CMACH or UZERO and the other will be computed within the 

code. 

IRSTIN - The station number to be read in for restart 

If IRSTIN = 0, not a restart case 

IRSTOT - The interval for saving restart information 

If IRSTOT = 0, no restart information saved 

JRSTIN - Logical file name from which restart information is to be read 

JRSTOT - Logical file name onto which restart information is to be 

written 

NFILE - The sequence number in JRSTIN of the desired restart infor- 

mation 

NSAVED - The number of restart stations saved on JRSTOT j must be 

initialized in inputs to the number of stations already 
written (and to be preserved) on JRSTOT. Nominally initialized 
NSAVED = 0. 

NOTE: By setting JRSTOT = JRSTIN and NEILE = NSAVED, one file can be used 

for both read and write information, without destroying the already saved 
information. 

YSL0T(2) - Initial nondimensionalized boundary layer thickness 

NS - Number of the last marching step to be taken 

NFIRST - Number of the first marching step whose geometry is to be 

specified (usually l) 

X(NFIRST) - Nondimensionalized location of first two marching steps 
X(NFIRST + l) 
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AP - Ratio of consecutive step sizes in the marching direction 

x(i) = x(i-l) + AP + [x(i-l) -x(i-2)] 


REY 

- Reynolds nTjmber based on UZERO and YZERO 

KTURB 

- 0 Laminar with no slip boimdary conditions 

- 1 Turbulent with wall function boundary conditions 

IPLOT 

- 1 Write plot file (Logical Unit 8) 

- 0 No plot file 

IC0EF(I,J) 

- Controls dump of terms being loaded into the block tri- 
diagonal matrix 


1=1 J=l,10 Momentum Eqs. 

1=2 J=l,3 Continuity Eq. 

IC0EF(1,20) 

- 1 Dumps the stored geometric data at start of equation 
writing subroutine coefg. 

COEFL: 


IGDMP 

- 2 Dtimps equations loaded into matrix inverter 

- 0 No dump 


There are three types of output generated by the PEPSIG code. First is 
the restart information that is written on a user specified data file. This 
information is used to permit a problem to be run as a series of short com- 
puter runs instead of rmning it all at once. The restart is exact in that 
eases run using the restart capability yield exactly the same results as 
those run all the way throiagh without using the restart. 

The second type of output is the printed output delivered on paper or 
microfiche. This output contains tables of parameters of interest at each 
calculation plane. A sample printed output has been anotated to explain its 
salient features and is included below. 

In three dimensional flow problems the flow characteristics can become 
very complicated. To make the results of this calculation procediire readily 
available to the user, some visualization techniques were required. Ber- 
nhard Anderson of NASA Lewis Research Center has headed up the development 
of a family of advanced graphics packages that can be readily coupled to two 
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and three dimensional fluid dynamics computer codes. The computer code 
developed in this investigation has been coupled into this NASA graphics 
system to automatically provide Calcomp, l6 mm film strip, and even movies 
of the results. This is the third form of output. These graphics capa- 
bilities can be used during code development to help highlight problem 
areas. The graphics are expecially useful in presenting the results of 
calculations to help the user understand how the features of his duct 
design affect the flow of gas within the duct. 
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tcpltmf test case raeph le^'^ 

r-T- CuRVFD CENTEPEINE 

Kr,TL^Pt.PSlS. scratch File 

te^EPSlG.ABSA 

m|ad:i|.m| 71‘*’ 

P( 2 ',l) = l‘. 0 » 

P{ 3 » 1 ) 50 * 0 ’ 

V 11 = 0 1 35 » 2 . 0 . 1 . u . 0 . . Stations 1-5 

|El|]|p.sTlN= 0 »If<^TOT= 5 rNFIEE 0 Write Restart File 
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V r 1 )= 0 » 0 » 
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CMACH=» 1 » 
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SEND 


Read Restart and 
Compute Stations 6-13 
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.31298+00 

.32143+00 

.31298*00 

■.29578 + 00 
.27010+00 
.25083+00 
.23701 + 00 

■.228B7 + 0D 
. 22619+00 


.27024+00 

^7T35B+"00" 

.28376*00 

:^ 8 i?St 88 

.35347+00 
.37872+00 
.38929+00 
.37872+00 
135347+00 ■ 
.32513+00 
.30109+00 
.28376+00 
.27358*00 

. 27024+00 


.31429 +00 
131530+00 
.33051+00 

;?li ?§:88 

.41417+00 

.44447+00 

.45714+00 

. 44447*00 

141417+00 

.38015+00 

.35125+00 

.33051+00 

'iwnm 


.3583^+00 .40238*00 .44643*00 


.37726+00 

.:SS^Ji:88 

.47486 + 00 
.51021 + 00 
. 52500+00 
.51021+00 
r474 8B*0D- 
.43518*00 
.40145+00 
.37726+00 


.42401+00 

,53555+00 
.57595+00 
.59286*00 
.57595+00 
15355S + 00 
.49020+00 
.45166+00 
.42401+00 
.4D773 + DD' 
.40238+00 


.47077+00 

.i$Hi :88 

-.59625+00 
.64170+00 
.66071+00 
.64170+00 
159025 + 00 
.54523+00 
.50187+00 
.47077*00 
-145245+00' 
.44643+00 


.lOOQQ.Ql 
.92857+00 
.85714 + 00 
178571+00 
.71429+00 
.64286+On 
.57143 + 00 
■ .50000+00 
.42657 + 00 
.35714+00 
.2857 1 + 00 
12T4 2V+00" 
.l4286+an 
.71429-01 
.DCOOO 


-i4iit3 

.49048+ 
.49716+ 
.51752+ 
155‘208 + 
.6L025+ 
.65695+ 

. 70744 + 

.72857+ 

.70744+ 

.65695+ 

. 60025 + 

155208* 

.51752+ 

. 49716 + 

.49048+ 


♦ OD"'—'. ^8 5T*00 ' 

00 .53452+00 

00 .54188 + 00 

00 .56427*00 

OD — 160228 + 00“ 
00 .65528*00 

00 +71764+00 

Do .77319+00 
00 ;79643*QD 

00 .77319+00 

DO .71764+00 
DO .65528+00 
00— '.•60228+00 “ 
DO .56427+00 
DO .54188+00 
00 .53452+00 


■ ■ "."Hrr+Tcr ' 

.57857+00 . 

.56660+00 . 

.61102+00 . 
165749+00 2 

:?^g?s:88 : 

.83893+00 . 

.86429*00 . 

.83893 + 00 . 

.77834*00 . 

.71030+00 . 

165749+00 — 2 
.6l}02 + 00 . 

.58660+00 . 

.57857+00 . 


.l285+TnT" 
62262+00 
63131+00 
6577?*00 
T02TCT+0D— 
76533+00 
83903+00 
90467*00 
93214+DD ‘ 
90467*00 
83903+00 
76533+00 
70270+00— 
65777*00 
63131+00 
62262*00 


— ;4Bddioi 

.66667+00 
.67603+00 
.70453 + 00 
1T5297+O0 ■ 
.82035*00 
.89973+00 
.97042+00 
.lDD0r*Cl 
.97042*00 
.89973*00 
.82035+00 
175291+W ■ 
.70453*00 
.67603+00 
.66667+00 



T£-sT-e*^- - R-« t Pn ttvif 7 =-Second Local 


r\ A T r w « < f-a » 

U M 1 1 U 1 XT3 1 D 


»**«- 


LEVEL 




ihet* 


Polar Coordinate 


24^ 




1 

2 

3 

4 

77lrr*DTr 

iitri+BD 

7 

0 

9 

10 


. . . 

T T 

«~D0 CO 

. 719 r=I!l 

■.“T979.0D 

■.716 37"mi 

.9 28^9017 

.SDDOiini — ;r7T99ira ■ 

~76979»mr 

lY 

Y 

.inoOC’Ol 18D.D0 

180.00 

leo.OQ 

180. QQ 

180.00 

180.00 

180.00 

180.00 1 

LSO.DO 

180.00 

"16T * 14 

16r« 1 4 

16 7*14 

16'/. 14 

16T.19 

■167;n 

~I6T.1-9 

B7.T9 ■ 

167,19 

13 

#8671*i^02 
.7857 1400 

159.29 

159.29 

159.29 

159.29 

159.29 

159.29 

159.29 

159.29 

59.29 

159.29 

191.93 

12 

141.43 

141.43 

141.43 

141.43 

141.43 

141.43 

191.93 

141.43 

41.43 

11 

.7m29^Q0 

128.57 

128.57 

128.57 

128.57 

128.57 

128.57 

128.57 

128.57 

26.57 

128.57 

lU 

i6‘IZ86»6^ " 

115.71 - 

115.71 

115.71 - 

115.71 

115.71 ~ 

115.71 

115.71 

115.71 — ■ 

15.71 - 

■ 115.71- 

9 

. 57143*00 

102.86 

102.66 

102.86 

102 .66 

102.86 

102.86 

102.86 

102.86 

02.86 

102. B6 

8 

•50000*00 

.42857*00 

9Q.QU 

90.00 

90.00 

90.00 

90.00 

90.00 

90.00 

90.00 

?9:SS 

90.00 

7 

77.19 

77.19 

77. 19 

77.19 

77.19 

77.19 

77.19 

77.19 

77.19 

8 

.357 IUffn 

6H .TT" 

69. 77 

69779 

~6;il .77 

— 69729' 


“69779 

69 .T9 

69729 

- 69729 - 

5 

.2857 1*00 

51.43 

51.43 

51.43 

51.43 

51.43 

51.43 

51 .43 

U:S? 

51.43 

51.43 

<1 

.ZltZ^^OO 

38.57 

38.57 

38.57 

38 .57 

38.57 

38.57 

38.57 

38.57 

38.57 

3 

r 

.m2e6«Do 

.71929-01 ' 

25.71 , 

■ 12.86' 

Uilk 

25.71 

12.86 

25.71 

iz;66 

-illlt ■ 


25.71 
12.86 

25.71 

12.86 

25.71 
12.86 " 

-U:U 

1 

•OnoCQ 

.00 

• 00 

• 00 

• 00 

• GO 

.00 

• 00 

.00 

• 00 

• 00 



IZ: 

11 

12 

13 

19 

.l§nr)ioi 


■7 T ' 

■■ .7I43*CC •' 

.7057+DTr 

.e5n*ao" 

.9286*00 

15 

14 

.100PQ‘01 

.9?857*D0 

IBC.CQ 

167.14 

160.00 

167.14 

leci.oo 

167,14 

18 0.00 
167.14 

180.00 

167.14 

13 

.85714.QM 

154,29 

154.29 

154,29 

159 .29 

154.29 

|i 


19 1.9 3“ 
IZ^.S! 

' “191.93 “ 
128.57 

191.93 " 

128.57 

“191.97 ■ 

128.57 

191.93 

128.57 

10 

9 

•64286*00 
.57 143*2'' 
.50orc+nn 

115,71 
lc2,86 
9". ''3 

9 0.00 

115,71 

102,86 

9G,0D 

115.71 

115.71 

- 

7 

.42857*00 

.35714*00 

77.19 

77.19 

77,14 

77.19 

77 . 19 

6 

64,79 

64,29 

64,29 

64.29 

64.29 

5 

.2857 l*nc 

51,43 

51.43 

51,43 

51 .43 

^.43 
JB757 

X 

■.21929*go 

-- 38.57 “ 

78.57 

“38.5 7 

- 38. 5T 

.19286100 

25.71 

25.71 

25.71 

25.71 

25.71 

2 

.71429-01 

12.66 

12.66 

1 2.86 

12.86 

12.86 

1 

•ocor: 

,no 

«0Q 

• no 

• 00 

• 00 


IjO 

00 
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LOCATION OF CENTERLl 


BWD = 

IDUCT = 1 


1-TH STATION 1 S- .OOrOODOQ 

.INE ( .aOOOOOOQ T .aOOOOOOQ < .00000000 I DX= • 

KOORD- 1 

IFRCT : 1 ITCNT r 2 VELOCITY BOUNDARY CONDITIONS 

. I H A C H - ^ vtuocilY '' vLCIOK SlSRL AYED IN CUMrS ^ A UonAl” ^ ^ t B 6R DI 8 A T ES 


oooooooo 

6 e 

6 t 


12 

12 


3 

20 


Initial Values of 


- X" 

level 1 ***** u-vEL Streamwise Velocity 




z : 

lY Y 

15 ;ixnnro*oi 

l*t . 92857*00 
13 .857 14 * 0 ; 

12 .7857 1*00 

~n "; 71429*00 
10 . 64286*00 

9 .57 14 3*00 

8 . 50000*00 

“7 mZflSTABO 

6 . 35714*00 

5 .2857 1 * 0 C 

4 . 21429 * 0 '’ 

" 3 " . 14285*00 
2 . 71429-01 

1 .00000 


1 

.0000 


ngoooYjr 
. 10000*01 
. iuOon*oi 
. 10000*01 
. 10000*01 ■■ 
. 10000*01 
. lunOL'*Ol 
. 10000*01 

■ 7}ooflc*nr‘ 
. icooo.ol 
. 10000*01 
. 10000*01 

■ .ioour*oi 
. 10000*11 
. 10000*01 


7 

. 7143-01 


. 1429*00 . 2143*00 . 2857*00 . 3571*00 . 4286*00 . 5000*00 


"V 

, 5714*00 


10 

. 6429*00 


^'TiraQ»vi 
. 10000*01 
. 10000*01 
.i''Lao*ai 
■.iru 03 ‘ 0 i 
. 10000*01 
. 116 JO*Ol 
. 10000*01 
10003 * 01 ' 
. 11000*01 
. 10000*01 
. 10000*01 
. 10500*01 
. l'’LU 0*01 
. 10600*01 


.10000*01 
. 10000*01 
. 10000*01 
. 11000+01 
. 10000*01 
. 10000*01 
. 1 0000*01 
. 10000*01 
'. 1 0000 * 01 ” 
. 10000*01 
. 1 0000*01 
. 10000*01 
. 10000*01 
. 10000*01 
. 10000*01 


.10000*01 . 1UQUD*U1 


. 10000*01 
. 10000*01 
. 10000*01 
iiooor*oi 
. 10000*01 
. 10000*01 
, iuunu +31 
■ 71 DCDD» 0 r 
. 10000*01 
. 10000*01 
. 10000*01 
. 10 CD 0*01 
, 10000*01 
. 10000*01 


, 10000*01 

. 10000*01 

‘.mm 

• ioQon*oi 
•looor+Di 
* 1000^*01 
-.iDOor^oi- 
. 10000*01 
•looor^oi 

•lUOOO^Ol 
. ioncp*ni 
• iuoor+01 
.looon^oi 




. 10000*01 

:mm 

. 10000*01 

. 10000*01 

. 10000*01 

■;iDooo*gr 

. 10000*01 

. 10000*01 

. 10000*01 

. 10000*01 

. 10000*01 

. 10000*01 


. 10000*01 
. 10000*01 
. 10000*01 
. 10000*01 
. 10000*81 
. 10000*01 
“;ioooo*oi 
. 10000*01 
. 10000*01 
. 10000*01 
. 10000*01 
. iuooo+01 
, 10000*01 


.49 6 ^ 3 *0 0 
, 99732*00 
. 99855*00 
. 99974*00 
. 10000 * 01 ■ 
. 10000*01 
. l 0 Q 00 * 0 l 
. 10000*01 
“J 10000 * 01 ' 
. 10000*01 
. 10000*01 
, 99974*00 
.' 99855*00 
. 99732*00 
. 99683*00 


. 97589*88 

. 98077*00 

. 98824*00 

.99514*00 

. 99993*00 

. 10000*01 

. 10000*01 

“: I9857 : 88 ' 

. 99614*00 

, 98824*00 

. 98077*00 

.97589*00 

. 97422*00 


.9169 7*1 


.91 


1986*00 
. 92873*00 
, 94395*00 
, 96497*00 

. 99994*00 

iimm 

• 96 M 97+00 

• 99395^00 

. 92873*00 

• 91986 « 0 Q 

• 91697«00 


— 

“ TZz 

T 1 

■ 12 ' 

13 

- m - 

T 5 ' 


2 : 

. 7143*00 

. 78 5 7 + CO 

. 8571 .00 

. 9286*00 

. 1000*01 

IS 

.inooc +01 

. 81588*00 

. 66779*00 

. 47552*00 

. 24797*00 

.OOOOr 

14 

. 92057 *Da 

.0 19 S 9 * no 

. 67161*00 

. 47869*00 

. 24976*00 

•CDOOO 

13 

i 857 I 4 ^' 0 P' 

"i 0 ll 2 D* 7 D 

V 68 37^*00 

, 48880*00 ■' 

•. 25553*00 * 

iroroc 

12 

. 78571*00 

. 05217 *nO 

. 70631*00 

. 51788*00 

. 26653*00 

.00000 

U 

. 71429*00 

. 88434 +nO 

. 74283*00 

. 53976*00 

. 28521*00 

.OOOOr 

10 

• 64206 +CO 

• 92642 ^T 0 

• 79 507 + OO 

• 5 R 875 *na 

, 31480*00 

•nuoor 

■ 9 

.57 14 3 +rjn 

• 96523 *’'Q 

. 65533*00 

. 5491 C ‘10 

. 35304+00 

.oocor 

6 

.soooc*or 

• 9 799 1 * 0 J 

. 80418 +OO 

. 68152*10 

. 37459*00 

.ouocr 

7 

.4285 7*00 

.9652 5 * 0 J 

• 85533*00 

• 6491 0*00 

• 353 U 4 +U 0 

.otocn 

6 

.357 l 4 > 0 r; 

• 92642 *CU 

. 7*3 587*00 

. 58875*00 

. 31480*00 

.OCOGP 

”5 

. 2857 1*00 

■ .B 8434 * no 

-.74 28 3*00 

.5 3976 * 10 ' 

. 28521 » 00 ~ 

.00000 ■■ 

4 

. 21429 +OG 

. 85217 +CO 

. 70631*00 

. 50708*00 

. 26653*00 

.ruGco 

3 

• l 42 P 6 *nti 

• 8312 L*nj 

. 60375*00 

• 46880*00 

• 25553 *un 

•ounon 

2 

. 71429-01 

.0195900 

. 67161*00 

. 47869*00 

. 24976*00 

.ocooo 

I 

.r'^onr 

. 0153^00 

. 66779 *Cn 

. 4 755 2*50 

.I 4797 *CT 

.OLor.i 


I 
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level 1 •*♦** 



IZ- 

1 

2 

3 

tY 

7 = 

Y 

.none 

.7143-01 

• 1429 

lb 

.loauc^Di 

.OOD'OO 

• cnoou 

■UOOOu 

IM 

.9e857*oa 

•oouou 

•incoo 

•OOQDO 

13 

• 85714+OC! 

• ?coj: 

• w 0 0 0 0 

.00000 

12 

.7657 1*00 

.'■.uroo 


.cua:io 

11 

. 7 1 4 9 ♦ G 1 

• *. L L' o 

• u ■'cOO 

.nnrer 

I'J 

•64286+00 

. OOCI JO 


• nenoo 

9 

.5714 3*0'? 

• ouuc: 

• UCOiJO 

,r,pqon 

8 

•snocc+co ■ 

• Conor' 

.001.00 

.nmun 

.oonoo 

7 

• 4285 7 + 00 

• cunoc 

■'ianono 

6 

.357 14*00 

.olcj'j 


.'70000 

5 

.28571+00 

• u i.' U U f.' 

. JO.jOO 

•nauoo 

4 

.2 14 29 + 00 

• CLDOU 

•enuoa 

.nono j 


. 14 236*0'3 

.Cor jLi 

. iJOOQj 

. JCOOC 

2 

.71429-01 

• cooor; 

•uToOO 

.niOOO 

1 

..inooo 

• naojc 

•bCwOU 

. 0 jnon 


12; 

Z ; 

.iaooo«ai 

.92857*nn 
.857H*0n 
.7B571 + aO 
.7l429-*Qr 
.6<I286*0C 
.t,7l‘i8‘00 
.50000*00 
.M285 7*00 
_,.357 m»0Q 


n 

.7143»DD 
.coccio 
.cooar: 
'.Ooooc ■' ■ 

.OuHOO 
.00000 
.ojL'ao 
.QoCoO ' 

:S8ggf! 


12 

. / 857*00 

« oncoo 
•uTlOO 
•jOuOO 
•unoou 

.aTfiOO 
• CQUOO 
.Ij^coo 


•cooou 

•coocn 

•DUOOQ 

.njDCu 


•UULOO 
•UOUOO 
• jniiOO 


Lev€l 1 ««*** 


uo 

VO 




I 


.. .pfpsie ewvtD 

C£NTE»tlNE 

-TtST 

tEVY--- • • • 


\z- 

1 

2 

3 

4 

1 

[r 

2 - 

Y 

• nooQ 

. 7m3-Oi 

• 1429+00 

•2143+00 

! 

iS 

nT?OffTJ^I 

♦92857+00 

TTiffjac 

•ruooo 

• Li3ljUU 

• oO UOO 

.naOdO 

.00000 

• OCOOCi 

1 

13 

.8571t+0T 

.OOCdO 

.00 000 

.00000 

.00000 

■i 

If 

.7857 1*00 
.7l*t79*3D 

•Mm - 

•mm ■ 

•QGOOO 

.nnnnQ 

•nonoo 

• OODGO 

• uOGOO 

: 

LU 

•642e6«00 

• UOOQQ 

• OCuGO 

•QOCOO 

9 

6 

• 57 !«♦ 3*00 
•snoDC^oo 

•ooour 
. ojoor 

•acLOo 

•onujo 

•°C00Q 
• nriGOO 

:8?g88 


.42857*00 

.35714*00 

>00000 

•oocoo 

.00000“ 

.00000 

6 

•CoOCU 

:gfl^88 

•OGOQO 

•cnnoo 

.00000 

5 

.2857 1*00 

• OLinOO 

.00000 


• 2m29*30 

■ COCJO 

■ UG (jO 0 

.00000 

•OOGGO 


T" 

;i42B5*DC 

• noHCiD 

-.dPLOC 

.00000 ■ 

■.00000 

2 

.71429-01 

.00000 

.coloo 

.00000 

.ouuoo 

1 

•conoo 

• ouau'] 

•OOLOJ 

.OOOOD 

•00000 




• IS 

,iono*Gi 

.color 

.ocour 

vocnocT 

.Donoo 

.ooQon 

'MW 

• OOOtJO 

•cuoon 
.ouooci 

'MW 

•OOOOP 

.00000 

•ronon 


Initial Values of 
Radial Velocity 
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ttAJ£_ftl_U4aA. 
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l-TH STATION , S: .QODDDCIDO 

LOCATION OF CENTERLINE ( .00000000 i .00000000 . .00000000 

K00R0= 1 

BWD : 0. IFRCT - 1 ITCNT = 2 

IDUC T : 1 IHACH_:_J, N OHACH r J 

I 


I 




ox: .00000000 


Level i ***•• rho •••** 


Initial Values of 

Density 



I2 = 

1 

2 

3 

4 

5 

6 

7 

« 

9 

10 



■ vooro ■ ■ 

■ .7183*0T 

— .'14 29 +110 

.ZIUS+TTO 

;7B*57 + on 

.357I*on 

.828t*00 

"'iSDOO^UU “ i5Tt4*inr~ •B479TTO 

lY 

15 

Y 

•lOQoa^oi 

•loooQ^ai 

• 10C'3Q«01 

•iaoao+01 

.10000*01 

• lODon+ai 

.10000*01 

•10000*01 

.99999*00 .99990*i)0 

•99968*00 

lA .’T2O5Tnj0r 
13 .«57m*on 

•100UU*0I 
• 10000*01 

:\mm^ 

:WTSr 

:188Sni 

."10000*01- 

.l000''*01 

. jooou*in ' 
.10000*01 


• 9V999 ?W9T0*Tin 7 *Qu 

.99999*00 .99992*00 .99973*00 

12 

.7857 1*00 

• lQPO'J + 01 

.10000*01 

•logoo+Qi 
. iQOoo*ni 

•1Q000«0} 

.10000*01 

.10000*01 

• 10000*01 
• Icooo*oI 

: 18888:81 

1 .99995*00 

.99978*00 

11 

#7m 29 + 00 

• XuOCiO^Ol 

• inuou+Dl 

• laooo+c! 

. ionor+oi 

•IODOO*Qi 

[ .99998+00 

.99986*00 

■ 10" 

•6U28ei*nn 

.TL00C*DI 

■'110000*01 

.ICDDD^Ol ■ 

.10000*01 

.lCOnr-*01 

.10000*01 

.10000*01 

.10000*01 

1 "".10000*01 

.99995*00 

9 

8 

.57IA3*0D 

.5CUUD*D0 

. 10000*01 
.iDour*oi 

. 10000*01 
.10000*01 

•icngo+gi 
.1 Ciooo*ai 

.10000*01 

.10000*01 

•lUOOP^Ql 

•luDon^oi 

:i8888:8} 

:}8888:8} 

: 18888:8] 

1 .10000*0) 

1 .10000*01 

•99999*00 

•10000*01 

7 

.‘(2857*D0 

.IQD'jri*Ql 

. lOCOO^Ol 


.10000*01 

• iooor+-oi 

.10000*01 

.10000*01 

.10000*01 

1 .10000*01 

•99999*00 

% 

.35717**00 

' .’10000*01 

• incoo*oi 

.1000D*DI' 

■.10000*01 ■ 

.10000*5} ■ 
. iooor*nI 

,10003*01 ■ 

;}oogo+Di 

.10000*01 

;1000D*03 

1 .UODDO*!!! 
1 .99998*00 

.77995*0(1 “ 
•99986*00 

5 

.2857 l^-OO 

. 10000*01 

. 10000*01 

.10000*01 

.10000*01 

•IOOOU*Q1 

.10000*0] 

•loooQ+oi 

H 

.21429+00 


. ioano*oi 

•inooo +01 

• lODQO + Gl 

• looon+Qi 

.iuoor*oi 

:i8888:8i 

• 10000*01 
.10000*01 

[ .99995*00 

.99978*00 

.99973*00 

3 

.18286*00 

•iOLOQ^UI 

.10000*01 

.10000*01 

•99999*0( 

3 .99992*00 

■ 7" 

.71829-01 

. 10000*01 

.10C0D*01 

.10000*01 

• I OuH 0 * G 1 

.1COCF*01 

.10000*0} 

.10000*01 

.99999*01 

J .99990*00 

.99969*00 

1 

•unooc 

. U3fjui.'*ai 

• imjoo^oi 

, 1 Df^OQ + Ol 

. 1UD0(’*U1 

.10000*01 

•1000U*Dl 

• lonoQ+oi 

•99999*0( 

J .99990*00 

.99968*00 



12 = 

11 

12 

13 

14 

•iSro+ui 


• 7 - 

■ '.Tlo^^Tin 

.7857»riO 

•99889+00 

.8571*00 

.9ZB6iD3 


•ioaoc*oi 

•99933+Cm 

.99846+OQ 

.99813*Q0 

.99801+00 

14 

.92857+nn 

.999 3 u+ PQ 

.99890*00 

•99046*00 

.99613*U0 

•99801*00 

13 

— rz " 

.857 14*nn 
■.7857-l*0C 

•99930*OJ 
■~;99945*0a ~ 

•99894+CO 
■•9990O*OtI' ' 

.99848+00 
.99852*00 ■ 



11 

.71429*0P 

.99956+Pj 

.99911+00 

,99859+nO 

.99817*00 

.99801*00 

ic 

.642e6*0r< 

.99972+00 

.99927*00 

.99070*00 

.99820*00 

•99801*00 

.99801+00 

9 

.57143+00 

•99986+io 

.99946+00 

.99885+00 

,99825+UO 

- B' 

.5COCC*nO 

.9999*^ + 10 

.99956*00 

,9989 3*Oi3 

.99829*00 

.99801+00 

7 

.4285 7+on 

.99906*'^^ 

• 99y4fa + 0C 

.99885*00 

. 99825*00 

.99801+00 

6 

.35714+00 

.99972+Oa 

.99927*00 

.99870*00 

.99820*00 

.99001+00 

5 

.2P57 1 + 0'“ 

.99956+ GO 

.99911+00 

•99859+00 

,99817+00 

.99801*00 


■ .214Z9*0P 

.99945*00 ■ 

.99900*00 

.99852*00 

'.79815*00' 

.99801*00 

3 

. 1420 6*00 

.99938*00 

.99094*00 

,99848*00 

•99814+00 

.99801*00 

2 

.71429-01 

.99939+nj 

•99890+00 

,99846*00 

.99813+00 

. 99801*00 

1 

. ij 0 u c n 

• 999 3 3 + PO 

.99689+00 

,99846+QO 

.99813+00 

•99801+nu 


■p- 

M 



UP iT 0X10*8 




, PCP^Tts - cuRvttj— crnTe(ftrNf--TtST"C*5t — »Atr*rt-cmf 


«p 


I-TH station f S: .OQtlOOOOO 

LOCATION OF CENTERLINE ( .00000000 , .00000000 . .00000000 . » OX: .00000000 

NOORD= I 

BND : 0. IFRCT - I ITCNT r 2 

.1D.UCT INA.CH, :_1 NPMACH_= 0 


Initial Values of 


LEVEL I ***** PRESS ***** 


Pressure 


iz = 

2 

13 .ssnu^oo 

}? 

10 .OMSbi+OO 
9 •D7143+0C 

8 •snaac*QO 

-i-'Mnim 

5 .2857 1»00 

H ,2l‘U9*0U 
3 .m2e6*00 

I ".TiiizR-ai 
I .t’fiO'JD 


I 

.0000 


■■ . Tin J-Dl 


.H2TF.D0 


. JINTFTJD ' 


•28FT»DD 


,114 5*4*02 .11454*02 

.TI*454*r5 .“71454 FOr' 

.71454*02 .71454*02 


.71454*02 
.7l4 54*-’2 
.71454*12 
. 71454*12 
.71454*02 

.71454+02 
•71454+T2 
.71454+H2 
. 7 145 4 + 15^ 
•7l454+?2 


•7l454+n2 

.71454+02 

.7l4«^4^G2' 
.7U54 + 02 
.71454*02 

.71454*02 
.7 1454*02 
.71454*02 
'. 71454*02' 
.71454*02 


.7 145 4*02 
".71 4? 4 ^ 
.71454*02 
.71454*32 
.71454*02 
.7 1454 *12 
.71454*02 
.71454*02 

.71454*02 
.71454*02 
.71454*02 
. 71454'*62 
.71454*02 


71454*02 

71454Fin" 

71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 
.71454*02 
.71454*62 
.71454*02 
. 71454*02 


.71454*02 

“.TllBTTFflZ 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

'.71454*02' 

.71454*02 


;357i*iro 
.71454*02 

■^TiTr54»Q2' 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

“ lii454*o| 

.71454*02 
.71454*02 
.71454*02 
'.71454*02 
.7 14 54* 02 


.4^ 


B6*O0 iSOODFOD .'SriTRinJ" 




Z7F01T" 


.7 145 4*02 

'"r71«54'FlR- 

.71454*02 

.71454 *02 

.71454*02 

.71454*02 " 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454 *02 


,71454*02 

: 8 ? 

.71454*02 ■ 
. 71454 *02 
. 71464*02 
. 71454*02 
“.'71454*02" 
,71454*02 
.71454*02 
. 71454*02 

:]lslJ:Sr 


.71454*02 

:?lSir* 8 ? 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

".'71454*0? 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 


45 4*02 

1454*81 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
,71454*02 

.71454*02 

.71454*02 

.71454*02 

'.71454*02 

.71454*02 


15 

14 

13 

TZ 

11 

10 

9 

B 

7 

6 

5 

■4 

3 

2 

1 


.1000C*0 1 
.92857*00 
.85 7 1 4*CC 
Z7B57I*00 
.71429*00 
.64286*00 
.57 14 3*00 

.5c3oc*oo 
.4285 7*00 
.35714*00 
.2857 1*C0 
".21427*00 
.14286*00 
.71429-01 

.onooc 


.7143*1111 

.71454*02 

.71454*02 

.71454*02 

.71454*DZ 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.*1454*02 

.71454*02 

.71454*02 

.71454* 02 
.71454*02 
.71454*02 


12 

'. 7857*00 
.71454*02 
.71454*02 
.71454*02 
171454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*02' ■ 
.71454*02 
. 71454*02 
.71454*02 


Aln'*arT 

.71454*02 

.71454*02 

.71454*02 

i7I454*TTZ ■ 

.7)454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.7)454*02" 

.71454*02 

.71454*02 

.71454*02 


14 

".92S6*gn 

.71454*02 

.7)454*02 

.71454*02 

“.'7145 4*02“ ■ 

.71454*02 

.71454*02 

. 71454*02 

.71454*0? 

.71454*02 

,71454*02 

.71454*02 

.71454*02' 

.71454*02 

. 71454*02 

.71454*02 


. .^0017*0 1 

.71454*02 

.71454*02 

.71454*02 

^71454*02 ' 

.71454*02 

.71454*02 

,71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 



»»£PSi6- CURVEB— ££NTE<mf<E--^E-5^- £*«- Rdtf^ttVT 


LEVEL 1 ****•, CP/2 


^ — g» TE 01107 8 — WrtC ET ^ 

Initial Values of 
Pressure Coefficient 


X2Z 
7 “"i 


1 

TOD on 


.7141^01 


3 

■.74Z’5'*nO 


4 5 _ 6 _ 7 » I t . ■ 

T2T4J+W .7B57*nr ' .357T*D0 •^ZII6*bl) •SOdO'^QO isrlV + do 



.loSoD«01 

V 

•64634-06 

.54634-05 

• 6463‘t = Od 

.64634-06 

.64634-06 

.6<t53‘*-05 


.9285 7»D0 


,5'4634-1J5 

■.'6a-6T^a6'“ 

T64634=Ur 

.64634-06^ 

"1555 PC -'06 

13 

.85714+00 

•64634-06 

.54534-05 

•6Q63M-06 

.64634-06 

.64634-06 

.64534-05 

12 

.7857 1*00 

•64634-C6 

.54534-05 

•6463M-06 

•64634-06 

.64634-06 

.54634-05 

11 

.71429*00 

•64634-06 

.54534-05 

• 6463‘t-06 

•64634-06 

•64634-06 

.54534-05 

"TO' 

“'154 28 5 = 00 

■•‘6t»634-n6 

'.54634-06 

• 6463A = i,6 

.64634-06 

.64634-06 

.54534-05 

. 9 

.57143*00 

•64634-36 

.54634-06 

. 6463*4-06 

.64634-06 

.64634-06 

.54534-05 

8 

,soooc*oo 

.64634-06 

.54634-05 

.64634-Q6 

,64634-06 

•64634-06 

.54634-06 

7 

.42857*00 

•64634-06 

.54534-05 

.64634-06 

•64634-06 

•64634-06 

.64534-05 

— r 

"".TSlI^mO — T57i 534=05“ 

-154 534=06" 

■.6467*r=XT5- 

~.6463‘*=06" 

".6 4634-06- 

"155534 = 05 

5 

.28571*00 

.54634-05 

.54b34-06 

.64634-06 

•64634-06 

•64634-06 

.54634-05 

4 

.21429*00 

.54534-06 

.64634-05 

•64634-06 

.64634-06 

.64634-06 

.64634-05 

_3 

.14285*00 

.54534-05 

.54634-0^ 

•64634-06 

.64634-06 

.64634-06 

.64534-06 


-“Th4Z9-fTI ' 

".54634-06 ' 

154634-06 

.64634=06 

.64634-06 

.64634-06 

.64634-06 

1 

.00000 

. 54534-05 

.54634-06 

.64634-06 

.64634-06 

.64634-06 

.64539-05 






-•30734-Oe* 

-•30734-06 

-.30734-06 


-.30734=06 
•64634-06 
•64634-06 
•64634-06 
iTOTJl 

:; 38 ? 3 J-- 8 g 

-.30734-06 

-.30734-06 

-.30734-06 


-.30734-06 

'^i3D734-06 

-.30734-06 

:; 58 l? 4 ^- 8 ^ 

4=06- — 

::J 8 ??!: 8 I 

-.30734-06 
-.50734=06 ■ 

-.30734-06 


VW-- 


15 .100DC*ai 

14 .V2857*or 

13 .8 5714.00 

iri7f4”2^:8^“ 

10 .64286.00 

9 .57 14 3.00 

“ B "'.50000.00 
7 .4285 7.00 

6 .357 14.00 

5 .2857 1.00 

2 .71429-01 

^1 . .OOODD 


11 _ 

■ iTiBjior 
.30734-06 - 
.30734-06 - 
_i_30 73«-06 - 
. J(r734 = ff'6 
. 30734-06 
30734-f 
.30734-06 
. 30734-06 
. 30734-06 
.30734-06 
. 30734-06 
”. 30734 = 06 
.30734-06 
. 30734-06 
. 30734-06 


12 
” -.78 
.3073 
.3073 
. 3073 
'.3073 

.3073 

.30734-06 -.3073 
.3073 
1 3073 


.3073 
;3073 
. 3073 
.3^73 
.3073 


STiffir 

‘♦-06 

4-06 

4-06 

4 = 06" 

4-06 

4-06 

4-06 

4-06 

4-06 

4-06 

4-06 

4-06“ 

4-06 

4-06 

4-06 


13 

" .u^n^nci 

.30734-06 ■ 
.30734-06 • 
.3C734-06 - 
.3 j7rir=n6" ■ 
.30734-06 - 
.30734-06 - 
.3 u734-T6 - 
.30734-06 ■ 
.3C734-06 - 
.30734-06 - 
.30734-06 - 
.3073*J=06" = 
.30734-06 • 
.30734-06 ■ 
.30734-06 - 


14 

.30734-06 - 
.3C734-06 - 
.30734-06 - 
•3c734=q6 • 
.30734-06 - 
.30734-06 • 
.30734-06 - 
.30734-06 ■ 
.30734-06 • 
30734-06 - 
30734-06 - 


A^r 


.3n734-g6 

.30734-0 


6 - 


.30734-06 

.30734-06 


Pff+Ol 
.30734-06 
.30734-06 
.30734-06 
.30734=U6"^ 
. 3Q734-C6 
.30734-06 
.3C734-n6 
.30734-56 
.30734-06 
. 30734-06 
.30734-06 
.307 34=756 
.30734-06 
.30734-06 
30734-06 


LEVEL 


HaCH »♦*** 


4=- 

LO 



-p- 

p' 


-R«»sie— cOTvee-e^NTEKt-me- -^ e sT-e*se- ee v-y— - 


IV 

IE 

TV 

13 

12 

11 

10“ 

9 

S 

7 


IZ = 

Z : 


.lOOQCtOl 

-;97rfTPiKr- 


1 

.oono 


.7143-01 


.1429*00 


.2143*00 


.2857*00 .3571*00 


.85714*00 
.7857 1*00 
.71429*00 
■.64286*00 
.57 14 3*00 
.5nonG*op 

^.4285 7*00 

5 .'35 7 r4*u'rt 

5 .2857 1*00 

4 .21429*00 

3 .14286*00 

T .7I429-TJ1 - 

1 .00000 


ICOOC* 00 .10,000*00 .13000*00 

lijoiin+flB TiO 007371)7 rroQ57'*7D 
10000*00 . 10000*00 . 10000*00 


. 10000*00 

:ii^88g:33 
. 10000*00 
.100CC*00 
.10000*00 .10000*00 

:}88SS^?n]8888:'’‘’ 


. lOOOC + QO 
•lODOO^QO 
;ionoc*f73 
• ioouu*no 


. lOOCO + 00 . , - . 

ninaou+nn .IQOOQ^QO 

•laaoo^oo .tnoan+no 

1?88§:88 :!8S88:S8 
.10000*00 
.1 OODUITJU 
. 10000*00 
. 10000*00 
.10000*00 
. 10000*00 
,10000*00 


. 10000*00 
.13000*00 
.13000*00 
. inooQ*oo 


.10000.00 tlOOOCDO .10003*00 

0 — TiTDO 0770” ;o ooon*Dtr 

0 . 10000*00 . 10000*00 

. 10000*00 . 10000*00 . 10000*00 

.10000*00 .lOQQC^OO 

‘. 10000*00 ‘ .ioocr*co 
. 10000*00 . 10000*00 
.1QU00*U0 .10003*00 

.10000*00 .10030*00 

”.'3000 0*00' ‘ ■ 

. 10000*00 

. 10000*00 . 10000*00 
. 10000*00 . 10000*00 
. 10000*00 . 10000 * 00 ' 
. 10000*00 . 10000*00 


•10000*00 
TT70ff0* DO 
. 10000*00 
. 10000*00 
. 10000*00 
. 10000*00 

:i8888:88 

.10000*00 


:}8888:8H :}8888:88 


•iDUOa^OQ 

•iDUDU*00 

•XDDOO+OD 

•lOUOO+OO 


.4ZSb«QD 

•lonoo ♦no 
7100DD»0a" 
•10000*00 
• 10000*00 

• iongo*oo 

• 10000*00 

:}8888:88 
• 10000*00 
■.10000*00“ 
. 10000*00 
.10000*00 
. 10000*00 
.10000*00 
. 10000*00 


-0«TE- 01107 a— 
8 


Pt a r- r- 

f" W UL 




•5000*00 


*-01 


. 9969 2-01 
“,W7317 
.99854- 
. 99974-01 
. 10000*00 
. 10000*00 

• )Q09Q*Q0 
. lOooo.oo 
. 10000*00 
V1D00D*DD 
. 10000*00 
.99974-01 
.99854-01 
.99731-01 
.99682-01 


.5714*00 


.97417-C 



10 

.6429*00 
.91 683-01 

„ — i^9Tr-[rr“ 
ox .92860-01 
.94385-01 
.96491-0 
.98584-0 
.99860-0 
. 99994-0 
.99860-01 

.9868^7 

,96491- 
.94385- 
.92860- 
.91972-01 
.91683-01 



12 = 

u 

12 

13 

14 

15 


- 7 = ■ 


■■ .■78E7*TTD 

.8571*00 

.928SiOn 

.loro* 

15 

.ioonc*oi 

•81561-ax 

.66 742-01 

.4751 5-01 

.24773-01 

.ooooF 

14 

,9265 7*00 

•819i2-ai 

.67124-01 

,47832-01 

.24953-01 

.00000 

13 

.85714*03 

• 8 3094-ni 

.68339-01 

.48842-01 

.25530-01 

.Gcnoo 

.00000” 


■“.“7857 1*00 


“■27D595*0r 

■;5g75i‘*in ■■ 
.53938-01 

■^266Z9-or” 

!i 

.71429*00 

.88415-01 

.74249-01 

.28495-01 

.00000 

ID 

.64286*03 

•926Z9-11 

.79558-01 

.58837-01 

.31452-01 

.00000 

9 

.57 14 3*0(1 

.96516-01 

.89510-01 

.64872-01 

.35273-01 

.00000 

8 

.50000*0’' 

.97987-11 

‘.88399-01 

.6811 5-01 

.37427-01 

.OLODP' 

7 

.4285 7*00 

.96516-01 

.85510-01 

,54872-01 

.35273-01 

.00000 

6 

.35714*00 

.92629-01 

.79558-01 

.58837-01 

. 3 1452-01 

.OCOQO 

5. 

.2857 1*00 

•80415-ai 

. 74249-01 

.5 3938-01 


.00000 

4 

—.21429*00 

-.93 19**- 01“ 

".70595-01 

■.50751*01 

.00000 • 

■ 3 

.14286*00 

•03094-01 

.68 339-01 

.48842-01 

.25530-01 

.00000 

2 

• 71429-0 1 

.81932-01 

.67124-01 

.47832-01 

•24953-Ul 

• oco on 

1 

.a(53QO_ _ 

.81561-01 

,65742-01 

.47516-01 

•24773-01 

• OOOCIO 


Initial Val'ues of 
Mach Number 


P-EP5I0 --fWVEO-CENTESLm -TEST -CitE tEVV 


0*TE -0«778 


-P* «- 




InTEGR*TEO properties *T ST4TI0N 2 


now ivrr^ 


.'STTP 50. FT ; 


LOCAL TO ORIGINAL AREA RATIO 
LOCAL TO ORIGINAL MASS FLUX RATIO 


— Pin ST TLinr= — 

i:8888 


Check on Mass Conservation 

;3r7V¥7jjuBH;7-5ECi ^ 



-PEPSie - tttfrvte eEfETEWL-tNE-^E5T--€*5E R*tPH-tEV¥ 


/ 


Results of 1st Marching Step 


n * ♦ r it< < w» ■ 


BUD r 0. 
IDUCT = 1 


2-TH STATION , S= .00000000 

LOCATION OF CENTERLINE i .OQOOOOnO . .00000000 ♦ .NOOOOOOO^ 00> DX= • ilOOOOOOO^OO 

KOORD= 1 

IFRCT r I ITCNT : 2 VELOCITY BOUNDARY CONDITIONS 6 6 12 

- - JFRCT = I DENSITY BOUNDARY CONDITIONS 6 ij 12 


IMACH 


DENSITY BOUNDARY C0NDITI0^ 


3 

20 


-PAft g - 




LEVEL 1 


RADIUS 


IT Y 

15 .ioo oc«qi 

'"1. i92B!i7*tn5 

13 .8571U*on 

.7e571«0D 

.7mz9«cio 

.57 14 3*00 
.5anoc»oo 
.42857.00 
.35714.00 

:i!5?r.88 

.14288.00 
.71429-01 
.OUOOC 


12 

11 

16 

9 


1 

VtIDOD ■ 



:i8§88:3i 


50GOC-01 

.5Q0JC-nt 

.50000-01 


.7143-01 


.94048-01 
'774 715^r 
.95752-01 
. 10U21»00 
.10503*00 
.11069.00 
.11574.00 
.11786*00 
.11574*00 
.11059*00 
.10503*00 
. 10021*00 
.96752-01 
.94 716-01 
.94048-01 


.14 29-i'DO' .Z143*nO 


.18214*00 
“.“I8415POU 
.19025*00 
.20052*00 
21508 *00 


.1361 0*00 
'.1 394 3*00 
.14350*00 
.15042*00 
.15005*00 

■Aiuim 

•ia571*QO 
, J8iM9*00 
.171394D0 ' 

.j4350*00 
.13943*00 _ _ 

.13810*00 .18214*00 


.25357*00 
.24723*00 
.23208*00 
.21508*. 
.20062*01 
,19026*00 
.18415*00 


.zS 


57*00 


.22519*00 
".ZZB8T*0D' 
.23701*00 
.25083*00 
.27010*00 
. 29278*00 
.31298*00 
.32143*00 
.31298*00 
.29278*00’ 

88 

.23701*00 
.22887.00 
.22519*00 


3 I 7 I* 


on 


.428 6*00" ■ ; 5000*710 V5T1T*00" 


.6429*00 


.27024*00 

"."27758*00" 

.28376*00 

.30104*00 

.32513*00 

,35347*00 

.37872*00 

.38929*00 

.37872*00 

.35347*00 

.32513*00 
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.72519-03 

.81699-03 

.85280-03 

“♦B167r-*03“' 

.72599-03 

.69968-03 

.59199-03 

.55993*03 


:;l3?il:8i : 


-.59213 “ “ 



-.555?3-g3 


03 

03 -.61657 . 

!D3*-.72539**T)3 
-.88591-03 
-.10576-02 
-. 1 1287-02 

3 -188692-0 
03 -172503-03 
03 -.61669-03 
03 ■=155619=03 
03 -.52632-03 
03 -.51659-03 


5=0T* 
92231-03 
-.95156-03 
-.51963-03 
“*.63956*1)3 
-.85307-03 
-.11060-02 
2181-02 



399-g 
*.63927-03 

^:U\Vai 

-.92225-03 

-.91270-03 


15 


2 : 

.10000*01 


19 .92857* 00 ■ 

— 13 — .“85719*tITr* 


12 

li 

9 

8 

7 

6 


“ .“857 ] 

.7857 1*00 
.7l929*fir ■ 
. 69286*00 '• 
.5719!*0r, •- 
.50000*00 • 
.9285 7*00 ■ 
.35714.QU • 
17857 I»C0* 


9 .21929*00 

i :)H?li:8? 

““ .50000 


11 

.7193*00 
-.27176-03 -. 
.27950-03 -. 
.30356-03“ -1 
.35779-03 -. 
.97367-03 -. 
.69936-03 -. 
-.10199-02 -. 
-.11753-02 -. 
-.10192-02 -. 
-.69973- 03 - . 
*187395*03“=^; 
-.35775-03 -. 
-.30376-03 -. 
-.27995-03 -. 
-.27169-03 *. 


17“ - 

.7857*00 
12355-03 
12879-03 
19 979*1)3 “ 
18 399-03 
27201-03 
96436-03 
77189-03 
99855-03 
77179-03 

18 396-03 - 
19513-03 - 
12871-03 - 
12 398-03 - 


.8^71*00 

-.77022-05 

■.10351-09 

•.r85rrr*i}ir 


— in 

.9286*00 
.19917- 


.TIoo *01 

0000 c 


-.39278-09 
-.91235-09 
-.21660-03 
*.99155-03 
-.58219-03 
.99 197-03 


09 

— *‘ 8 {j 88 & 

.00000 

:88888 
.OCDOP 
.OOOOC 

.ooooc 


.10953-09 
-.60575-05 
-. 99687-09 
-113958-03“ 
-. 19099-03 
-.13955-03 


X679-Q3 

-•M97H9-0M 

•QbOOr 





1 195*09“ 
932 3-09 
8581-09 
039 5-09 
7511-05 

-.60501-05 

.10930-09 

.16526-09 

.18792-09 

".19925*09 

.DOOOt! 

.00000 

.OOOOC 

.60000 

.00000 
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PCPSI 6"--€ WVEe - et-ftTefrtfWE TEST - eO - St — ttv y — - 


LOCATION OF CENTERLINE I .OODOaono 
BWD □. 

IDUCT_: .1 


2-th STATION t i- .OOCOOOCD 
’ OQOOOUOO 


KOORO: I 
IFRCT = 1 
IMACH r 1 


.ODOOOOOO'OO) 




DXr .AOODDOOO^OO 


-P*Gt- 


-ATfl- 


ITCNT - 2 

NOHACH - 0 


LEVEL 1 ***** RHO ***** 


12 = 
7 T 


lY V 

15 . 10000 »tn 

13 .85714»np 

12 •7857 1*00 

11 .71429*00 

In'" ;642B6*00 
9 .57 14 3*00 

8 .50000*00 

7 .4285 7*00 

■“■ S ' .■35 7 J 4* DT ' 

5 .2857 1*00 

4 .21429*00 

3 .14286*00 

2'“. 71429-01 

1 .00000 


1 

-.0000 ■ 


.10000*01 
TimmuTOT - 
. 10000*01 
. 10000*01 
.lCODO*Oi 
.99999*06 

,99999* 00 
.99998* ig 
.99999*00 
■V999 9 9+0CT' 
. ituuo*oi 
. 10000*01 
. 10000*01 
.iooor*oi 

. IDOuO* 01 


■ ; 7 i 4 T-tn 


.■n9 27»00 '^214T»no - ite^^TrtVr i357t*D0‘ ;42r6*Otr " ;sflo0*tra 2 S?lir* 1 Jl) IS^VrOtJ-- 


. 10000*01 
■ nTTOODAOT - 
.10000*01 
.10000*01 
.10000*01 
199999*00 
.99999*00 
.99999*00 
.99999*00 
i99999'*0ir 
.10000*01 
.10000*01 
.10000*01 
. idoo3*di' 
. 10000*01 


. 1000 0*01 . 1000 0*01 

•~1 UOOO ^ Ol - • TD • 

.innoo+ai •iaooo*oi 

— 10+01 •lOQQO+ai 

•iocao+01 

;99999* C 0 ■ 
•99999+00 
• 99999 +00 
•99999+00 
■'i99999*0Cr 
4 inooo+ai 
•loaoQ+ui 


•lOOOC 
•10000+01 
'•99999 + 10 

•99999+OU 

•99999+00 
•99999+10 
■;9979 9^+ ar ' 
•locoo+ol 
.10000+01 


•iDogg+gi .loooo+oi 
.iaoon+3f ■■ ; loooo+oi 
• irjpoa+ni •inoQo+ui 


ICOOO+Ql 
mOTTO^l}! ■ 
XGOOO+01 
lUQOO +01 
10000+01 
99999+00 
99999+00 
99999+00 
99999+00 
9 99 9? ♦'DO 
ioaor +01 
lOOOP+01 

iooor+01 

lODor+oi 

iuoon+01 


+ 10DD 0+01 
♦ IPODD+'DT* 
•lOQOQ+gi 
•loono+oi 
•10000+01 
.10010+01 
. 99999 + 00 
.99999+00 
•99999+00 
10000*01' 
0000+01 
+ 01 


1 ! 


0000 * 


. 10000*01 

;iD0D0*01 

.10000*01 


.10000*01 
'TiaODOFlTI 
.10000*01 
. 10000*01 
.10000*01 
.10000*01 
.49999 *00 
,99999*00 
.999«9 *00 

■■; ioooD*or 
. 10000*01 
. 10000*01 
. 10000*01 
. 10000*01 
. 10000*01 


.99999*00 .99990*00 .99968*00 


.99999*00 
. 10000*01 
.10000*01 
. 10000*01 
. 99999.00 

:II?9’?:88 

.TOOOD*Cri' 
. 10000*01 
. 10000*01 
.99999*00 
. 99999*00 
. 99999*00 


,99992*00 

.99995*00 

,99998+00 

;10000*01- 

. 10000*01 

,99999*00 

. 10000*01 

VlDODOFOl' 

,99998*00 

,99995*00 

,99992*00 

.99990*00 

.99990*00 


,99972+C 
,99978*t 
.99986*00 
-,9999i»»D0 ' 
.99999*00 
.1000 0*01 
,99999*00 
' i '9999 q *'0 ir “ 
.99986*00 
.99978*00 
.99972*00 
.99969*00 
.99968*00 



IZr 

2 - ' 

11 

• 7mr+DT3 

12 

■■ •785T+Tia 

. Hti +on 


15 

■ Tiapc+ui 

15 

.iOUOC+Ol 

•99933+^0 

•99689+nn 

•99890+0“ 

.99845+00 

.99813+00 

.99801+00 

14 

.92857+00 

•99934+nO 

.99846+00 

.99813+00 

.99801*00 

13 

.857 14+OC 

•99938+00 

•99893+00 

. 9 9 8i| 8 + 0 0 
• 998%2iaa'” 

.99814+00 

•99801+00 

* 17 

~.78 5 7r*or 

•99945+03 ■ 

.99956+70 

i999ao+ail “■ 

~.998l5+0Cr 

%998()l*0l) 

.99803+00 

11 

•71429+00 

.99910+00 

.99858+00 

.99017+00 

in 

*64286+00 

.099 7 1 ♦ 10 

.99926+00 

.99869+00 

.99820+00 

. 99801 + 00 

9 

.5714 j + on 

.9998 5+ no 

• 9994b+rj0 

.99884+00 

.99825+00 

.99801+00 

■ 8 ' 

■.5ncc+on 

.9999 1+ rO 

• 99955 + T3 

.99892*::0 

•99B2B+D1 

.99801+00 

7 

.42357+30 

• 999fi6 + '^iJ 

.99971+1J 

.99946+00 

.99804+00 

.99825+00 

•99820+QQ 

.99801+00 

6 

•35714+nu 

•99926+nO 

,99869+au 

.99801+00 

5 

• 28 5 7 1 + C'J 

• 99956+''U 

•9991 1+00 

.99858+00 

♦99817+UO 

.99803+00 

--If - 

".21429+00 

.99945+00 

.99900+00 

.99852+03 

.99848+00 

■'.99815+00'" 

'.99801+0C 

3 

.142P6*<30 

•. 09930*00 

•9969 3 *00 

.99814+00 

•99B01+0C 

2 

•71429-01 

.99934+ " J 

•9969U+P0 

.99846+10 

.99813+00 

.99801+00 

.99801+00 

1 

• LOUOC 

•99933+-0 

• 9<58fl? + UU 

, 998-4 5 + 03 

.99813+01 





, PTPSIG -CWRVEO trrrtRLlKC-TtJT C*St— R»tAtrt-EVY- -- 

2-TH STATION , Sr .aOPOOODO 

LOCATION OF CENTERLINE ( .00000000 t .00000000 , .UOOOOODO»00) 

KOOROr 1 

BWD r 0. IFRCT r 1 ITCNI r 2 

lOUCT : 1 IHACH : J NOHACH : 0 


0*ff-irtttT8- 

DX: .DOOOOOOOaOO 




LEVEL 1 


PRESS **•»* 


IZ = 
7 = 


1 

.0000 


. :iR3-01 


.IvZ^iDD .ZI#3*on .28F7iCP 


.3571t00 


.azs6«an 


.5D0iJ*13D "i 5715*00“ 


la 


IT r 

15 .10000*01 

K .“92 8 5 7*071“ 
13 *85714*00 

.71453*02 

•71453+02 

.71454*02 

•71454+U2 

,71454*02 

.71454*02 

-.“7745T.B2- 
.7145 3*02 

TTl4S3*02“ 

.71453*02 

.T145T»T]2- 

.71453*02 

— 71454W- 
.71454+02 

-.7TT57t*lI2 

.71454*02 

■.71454*02 

.71454*02 

12 

.7857 1*00 

.71453*02 

.71453*02 

.7 1453*02 

•71453+02 

.71453*02 

•71454+02 

iS 



.71453*02 

.71453*02 

.71453*02 

.71453*02 

.71453+02 

•71453+C2 

.71453*02 

.71453*02 

.71453+02 

.71453+02 

9 

.57143*00 

.71453*02 

.71453*02 

.71453*02 

.71453+02 

.7 1453*02 
.71453*02 

’Miim 

8 

.50000*00 

.71452*02 

.71452*02 

.71452*02 

.71452*02 

7 

. 42857 + 00 

.7i45:*a2 

.71453*02 

.71453*02 

.71453*02 

.7 1453*02 

.71453+02 

Tb' 

5 

■.357U+6g‘ 

.2857 1*00 

-'.7145 3*32“ 
.71453*32 

;7l453+02‘ 

.71453*02 

.71453*32 

.71453*02 

“■.71453*Q2‘ 

.71453*02 

:Ui!^5:8l 


i| 

.21429+00 

,7145 3*02 

.71453*02 

.71453*02 

.71453+02 

.71453*02 

.71454*02 

3 

•142B6+C0 

.71453*02 

.71453+02 

.71453*02 

.71453*02 

.71453*02 

.71454+02 

.71454+02 

.71454*02 


■ ‘.71459-01 

““.-71453*12 

;71453+0Z 

.71454+02 

.7 1454*02 

1 

•UDGOO 

. 7145 3+02 

.71453*02 

.71454*02 

• 71454^(I2 

.71454+02 

.71454*02 


,71554*02 

;7m54*D2"“ 

.71454*02 

,71454*02 

.71453*02 

.71453*02 

.71453*02 

.71453*02 

.71453*02 

.71454*02 

.71454*02 

.71454*02 

,71454*02 


.71454*02 
-.-71454*07 
. 71454 *02 
.71454*02 
.71454*02 
.71453*02 
.71453*02 
.71453*02 
. 71453*02 
.71453*02 
,71454*02 
.71454*02 
. 71454*02 

,71454*02 

.71454*02 


.71454*02 
-.“71454*112- 
,71454*02 
.71454*02 

:?}l!ir*8? 

.71453*02 
.71453*02 
.71453*02 ,71453*02 


.71454*07 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 


.71454*02 

-1-77454405 

.71454*02 

.71454*02 

I?i!l4*8l — 

.71453*02 

.71453*02 

71453*02 

mvs 

.71454*02 

.71454*02 


\l 

\h 

9 

-8 

7 

6 

5 

4“ 

3 

2 

1 


IZ: 

:m?^*8^ 

.85714*00 

“;7B571*00 

,71429*nn 

.64286*00 
.57 14 3*00 
.50.100*00 
.4285 7*00 
.357 14*1,10 
,2857 1*00 
-.“214 29*00 
.14286*00 
.71429-01 
.00000 


11 

.7143*00 


,7 145«*02 
.11464*02 
.71454*02 
■;TI4 54*02- 
.71454*12 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*12 
.71454*02 
,.71454*02 

.71454*02 
.71454* 02 
.71454*02 


12 

. 7857*00 
,71454*02 
,71454*02 
.71454*02 
771 45 4* or- 
.71454*02 
.71454*02 
,71454*02 
.71454*02 
.71454*02 
.71454*02 
.71464*02 
771454*07- 
.71454*02 
.71454*02 
.71454*02 


13 

' .857TF00 
.7 1454*02 
.71454*02 
.71454*02 
.71458*02““ 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*02 
.71454*02 - 
.71454*02 
*71454*02 
.71454*02 




2Sb*m 
. 71454*02 
. 71454*02 

. 71454*02 
. 71454*02 
. 71454*02 
.71454 *02 
. 71454*02 
. 71454*02 

. 71454*02 , . _ . 

“ 271454*02 “. 71454*02 

*71454*02 .71454*02 

.71454*02 -“ 

. 7I454*02 


IS 

*ioro*“Di 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 

.71454*02 


.71*54*02 

.71454*02 


v/1 


^•-pep^ig tuRVEtr sT-tfrst pttPt* ttvr- 


LEVEL 1 ***** CP/2 ***** 


■.7HT=T31 


.1429700 ^2ini3V0!r V255T7^Cr ,3tTl.DD .U2iTE7DIJ 25Dmjnnj ;5TTT»TJtr 2i?297Tnr 




.85714+aC 
.785 7 1*00 
.71429*00 
.64206*00 
.57 14 3*00 
.50000*00 
.4285 7*00 
‘.35714*00 
.2657 1*00 
.21429*00 
.14286*00 
.71429-01 
.00000 


-.47333-13 
-.53627-03 
-.66216-03 
-.00344-73 
-. 12925-02 
-.15357-02 
-.12944-02 

7.-70344-03 

-.66216-03 

-.53436-03 

-.47428-03 

7.44053-03 

-.43900-03 


• H7lj 147-03 

• 5 3 31U-03 
•65930-03 
.89295-03 
.,i?b87-02 
•15C33-02 

jill9970§" 

.65930-03 

.53246-03 

.47238-03 

.44281-03 

.43328-03 


-.42088-03 

•7.'4 77!rp03^ 

-.45330-03 

-.5 1529-03 

-.63736-03 

7.06910-03 

-.12429-02 

-.14775-02 

-.12439-02 

7.S6SlT70r 

-.63832-03 

-.51434-03 

-.45426-03 

-.42565-03 

-.41706-03 


-.39417-03 -.35794-03 -.31311-03 - 
7-,7finta7fl3”-. 36<r5T703‘‘^.32O7Ti-trr^ 


-.42469-03 - 
-.4847 7-03 - 
-.60208-03 - 
-.02810-03' - 
-.11953-02 - 
-.14209-02 - 
-ill943-02 - 

“262714 703 ■- 
-.60208-03 - 
-.48382-03 - 
-.42660-03 - 
7.39990-03 - 
-.39131-03 - 


-.38750-03 -.34268-03 
-.44377-03 -.39322-03 
-.55439-03 -.49622-03 
-.77103703 7,70221-03 


-.11276-02 -.10389-02 
-.13555-02 -.12592-02 
-.11276-02 -.10398-02 
-.77D8F7D3‘‘7.TDZ217D3 
-.55535-03 -.49717-03 
-.44281-03 -.39227-03 
-.38845-03 -.34268-03 
-. 36366-03 -.3 1788-03 
-.35507-03 -.31025-03 


-.25971-03 

-i2663B^“ 

-.28641-03 

-,334n9-03 

-.42851-03 

-.61924-03 

-.93205-03 

-.11390-02 

-.93300-03 

7".6173J-03 

-.42851-03 

-.33314-03 

-.28641-03 


-.22347-03 -.15671-03 
-.26638-03 -.19295-03 
-.35031-03 -.26352-03 
-.52101-03 7. 41039-03- 
-.80521-03 -.66025-03 
-.99404-03 -.82524-03 
-.80616-0 3 -. 65930-03 
7^52 101 ^ 3“=^Fn9937Tn" 
-.35126-03 -.26448-03 
-.26543-03 -.19104-03 
-.22347-03 -.15671-03 
7.20344703 -.14050703 
-.19676-03 -.13573-03 


-.93767-04 

-.12047-03 

-.17578-03 

=.29213703 

-.49431-03 

-.63355-03 

-.49526-03 

v;Z9nB=03- 

-.17578-03 

7. 8 137 0-04 ■ 
-.83277-04 


11 

■ .714J*iK[ 
-.44176-04 - 
-.38454-C4 - 
-.45131-04 - 
='. 6229 6 -54 = 
-.98536-04 - 
-.17674-13 - 
-.32170-03 - 
-.42374-13 - 
-.32265-13 ■ 
-.17674-13 - 
-.98536-14 ■ 
7.61342-04 = 
-.4513r-04 ■ 
-.38454-14 ■ 
-.4bl3r-"4 ■ 


12 

.78 5 T»-Ca 

.20334-04 
. 12705-04 
.15566-04 
-.24 14 9 704 ■ 
.43223-04 
.85 16 it-04 
. 16720-03 
'.22728-03 
. 16 720-03 
.85 184-04 
.43223-04 
.24 149-04 
. 15566-04 
.12705-04 
.2f334-04 


■ ".l5 7I*Tm 

.88904-05 

.31684-05 

.4 1220-15 

.5'9031-OS' 

.14612-04 

.29671-14 

.6 1342-34 

.84231-04 

.61342-14 

.29071-04 

.1461 5-14 

.69831705 

.41220-05 

.3 168 4-05 

.88904-05 


-.50757-05 - 
-.31684-05 - 
-.31684-05 - 
7-.4tZ2 07a5 ^ 
-.60294-05 - 
-.79367-05 - 
-.12705-04 - 
-.16520-04 - 
-.1Z7J5-04 - 
-. 79367-05 - 
-.61294-05 - 
7.31684-05 = 
-.31684-05 - 
-.31684-05 - 
-.4 122 0-05 - 


-.31684-05 

-.31684-05 

-.31684-05 

-.31684-05 

-.31684-05 

=.31664705 

-.31684-05 

-.31684-05 

-.31684-05 



LEVEL 1 ***** M»CH ***** 



.7H3*01 


.21*t3*00 .2857*00 


X3 .857Xi**00 
X2 .78571*00 
11 .71829*00 

■nT"”;6878 6^5TI' 
9 .57183*00 

? :5?|ra 

6 r35 718*00 

5 .2857 1*00 

8 .21829*00 

3 .18286*00 

— 2 rT18 29-m - 

. 1 .00000 


.loois^oo 

.10015*00 

H88i^:38- 
.10016*00 
.10015*00 
.10016 *00 
. iu.0 1 6* Ou 
.10016*00 
.10015*00 
.10015*00 
mroi5*w 
.10015*00 


.1001 8*00 

>inol5*oo !iooi8*ao -.ical8*oo .iooi8*oo 
■10015*00 .10018*00 .10018*00 .10013*00 

|10015*QD .10018*00 .10018*00 .10013*00 

iIoai5*OD~“.jDgj5*Dcr '.laontoo .iooi8*oo 
■10016*00 .10016*00 .10016*00 .10016*00 
•10015*00 .10016*00 .10016*00 .10016*00 

. 10016 *00 .10016* 00 ?10016*0g 510016*00 

no 01 6TDO .10015*00 — noo DO — n 00 it»dd 

.10015*00 .Iggi**oo . 10014*00 . 10013*00 

■10015*00 .10018*00 .10018*00 .10013*00 

■iocis*oo >iggi8*og .ioai8*oa .iooi8*oo 
rioal5*olr ■.looiTfBO .Ioqi‘i*do .iooi«*oo 
>10015*00 .10018*00 .10018*00 .10016*00 


.10018*00 .10018*00 


-. ical8*on 
.10018*00 


.10018*QQ 

■.iagi8*o9 

.10016*00 

.lgol6*gg 

.10016*00 


:mm^ 


.10013*00 

.10018*00 

.10016*00 

.10016*00 

.10016*00 


12: 11 12 

T '.TITUITTO ,7BOr*00 

15 .10000*01 .81852-01 .66618-gl .8 

18 .92857*00 .81881-01 .67L25-0I .8 

13 .85718.00 .83006-01 .68286-01 .8 

— ^.7B57T*00 — 1J5T5t- 0T — ^70501=01” .5 
11 .71829*00 .88313-Cl .78161-01 .5 

10 .682#6*00 .92520-01 .79i4lJ-01 .5 

’ •5Ii‘*3*Qg .96636-01 .85327-01 .6 

.5c9d6*00“". 9798 2-01 .88206-111" 16 

7 .82857*00 .96836-01 .85327-01 .6 

6 .35718*00 .92520-01 .79613-01 .5 

5 .2857 1*00 .88313-01 .78 181-01 .5 

•f ;2i4T9*oo"Tff5^io2-in — mreffi-in — .5 

3 .18286*00 .83006-01 .68286-01. .8 

2 .71829-01 .81881-01 .67025-01 .8 

i.._»0CigQQ„ _..818SJ-31__.65618-01 .8 


."Hti *00" 

7828-01 

7768-01 

6788 -01 

0696-01 

3868-01 

8709-01 

8682-01 

7820-01" 

6682-01 

8709-01 

3868-01 


8788-01 

7768-01 

7828-01 


18 15 - 

. V20 6 * 00 . luOO *0 1 

28779-01 .00000 

28976-01 .00000 

25566-01 .00000 


. 00000 “ 

.00000 

.00000 

.00 000 

.00000“ 

.00000 

.00000 

.00000 


.28535-01 . 

.31423-01 . 

.35088-01 . 

i37l«l-01 ■ . 
,35088*01 . 

.31623-01 . 

.28535-01 . 

*n6677-or~. 

.25566-0} , 

.28976-01 . 

.28779-01 . 


vn 

IjO 


I 


6 

3571*00 


7 

8286*00 


5000*00 


9 

5718*00 


6829*00 


.10018* 00 ,10011*00 

.10013* 00 T1001T*Og“ 
.iggi3*gg .10012*00 
.Io0l3*o0 .10012*00 


,99787-;0l .97815-01 


.10013*00 
.10014*D0 
. 10016*00 
.10016*00 
.10016*00 
n 00 19*00“ 


.10013*00 

.10018*00' 

.10016*00 

.10016*00 

.10016*00 


•99918-01 
.10005*00 
.10010*00 
. '10013*00 
.10015*00 

:i8si5^:as 


.98 065 - 01 
■98817-01 

“'110007*00 

.10018*00 


10019*00 — 110018*00— n 001 3950 “ 

i88iK88 :188ii:88 :0:gg 


»10013«QD 

»XDD13«OD 

aoom^oo 


•1QD12«QQ •99914-01 

iiOoiz*oo ' ^^9 789 -or 

►lOOll^OO 


^^9 789-01" 
•9974T-01 


“l985sl*oi 

.99913-01 


:1881ii:88 :iS§ir-g' 

TlDOOTTgo — 170650=0 
.99629.-01 .96831-0 

.98817-01 .98328-0 


.98065-01 .92803-01 

“;9T578 = rt — i9l7l5 = DJ 
.97815-01 


.91619-01 
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Check on Mass Conservation 
at Station 13 


PEPSI6 -tURVEO-CENTERLI-Ve-TtST— C*St R- * 1P H - 


tHT E- OnOT a H>G E Hfr 


INTE6RATE0 PROPERTIES »T STATION 13 


rx5w iurrs' ;izia6 sit.tt^ ;^j48i7-DiLBn./sti:i 

LOCAL TO ORIGINAL AREA RATIO liOOQR 

.LOCAL TO_.ORIGIN(|L M^SS FLUX RaTjO X.OOll 



Results of Calculation at 
Station 13 


PEPSIB CUP VEB-CE>tTE RUNE TF ST~c»5t- fr *L PH LEV Y- - B*TE- E tt tft T B E»*« 

13-TH STJTION , S- .it00168D5*Jl 

LOCATICN OF CENTERLINE ( .00000000 , .68 362240-01. .40000000* 01 1 DX r .50000000*00 

KOORD= 1 

BUD = 0. IFRCT - 1 ITCNT :: 0 VELOCITY BOUNDARY CONDITIONS 6 6 12 3 

IDUCT = 1 IMACH : 1 . JFRCT_3 1 DENSITY BOUNDARY CONDITIONS 6 11 12 2 0 


level 1 4»**» RADIUS ♦*»** 


IZ: 

- • Z ^ ' ' 

lY Y 
15 .10000*01 

T4 — rvrBSTTcnr 
13 .85714*00 

12 .78571*00 

11 .71424*30 

10 .64286*50 

9 .5714 3*00 

8 .5noco*or 
7 .4285 7 + 30 

E~~.35714*30 
5 .2857 1*00 

4 .21429*00 

3 .14286*00 

I 


1 

ttrotnr 



JOQD-ni' 
.5U0UC-01 

.icooc-01 

.50000-01 


“.7143*T7r- .14 79*00 

.94 04 8-01 .13 810* 00 

794 716 -nt .1 3W7Y-DO ' 

.96752-01 .14350*00 

.15042*00 
.16005*00 
■.17139.-1C 
.18l49*ij0 
.18571*00 

ni069 * 00-' + 

.16005*00 
.15042*00 
.14350*00 


; 2 ? 


4ITTT) '.ZBOTFOD .3571*00 .4286*00 


T5000*oti ;STi4*oo - -- r4r 




.13943*00 

.13810*00 


.18214 *00 

TiBTnBTtttr- 

.19026*00 

.20062*00 

.21508+UO 

:mTo 

.25357*00 

s24723*an 

723208*00- 

.21508*00 

.20062*00 

.19026*00 


.2261 9*00 
.72887*00 ■ 
.23701*00 
.25083*00 
.27010*00 
.2927P*0< 
.31298*U0 
. 32143*00 
. 31298*00 

-.29278*00 

.27Qir*C0 

.25083*00 

.23701*00 

.22887*00 

.22619*00 


27024*00 

2T3S8*00 ■ 

28376*00 

30104*00 

32513*00 

35347*00 

37872*00 

38929*00 

37872*00 

35347*00 

32513*00 

30104*00 

28376*00 

27358*00 

27024*00 


.31429*00 
-781 8 30*00 ■" 
.33051*00 

.iEA?i:88 

•41417+OD 

,M^^^7♦0D 

.MS7m^00 

,44*447^00 

• m4l7*D0 " 

»3601S^00 

.35l25^0Q 

t33D5l+OD 

*31B30^QD • 

.31429^00 


.36301 *00 . 4077 3*00 
.37726*00 .42401*00 
.40145*00 .45166*00 
.43518*00 .49020*00 
.47486*00 ■ . 53556*00 
.51021*00 .57595*00 
.52500*00 .59286*00 
.51C21*aB .57595*gg 


7474 86*0 


753556 *t 


.43518*00 .49020*00 
.40145*00 .45166*00 
.37726*00 .42401*00 

:3ll8i:88 


4-29*- tnj- 

•44643*00 
-745245*00 — 
.47077*00 
.50187*00 
.54523*00 
.59625*00 
.64170*00 
.66071*00 
.64170.00 
■759625*00 " 

:85iir*88 

.i*7077^D0 
•452B5»00 ■ 
.HH643^00 




a 

12 

13 

14 


■ ’ 2 r 

.T1»T*C0' 

■.7B57«0tr 

.S5 7i»an 

-. 92 8 6 + m 

15 

.10000*01 

.49048*00 

.53452*00 

.57857+00 

•62262*00 

14 

.92857*00 

.49716* 10 

•541B6+00 

•58660+au 

.63131*00 

13 

.85714*00 

.51752*00 

.56427*00 

.61102*00 

.65777*00 


■■■.T8571*tJU 

. 55208*00 ■" 

‘;6022B*00'’ 

-.6 5 24-9* aO ■ 

70Z7D'*QCr* 

il 

•71429+on 

•60025*10 

.65528+00 

.71030*00 

• 7653 3*00 

10 

.64286*00 

.65695*10 

.71764*00 

.77834*00 

,03903*00 


.5714 

•B^nnc+nr 

.428 5 7*'jti 

.3571M*no 
.2857 I + CJT 
■".21429 + 00 
.14286+OC 
.71429-01 
.uPfJGC 


.70744+00 

.72857*00 

.7U74 4+n,j 

.656V5*;ia 

.6U025+i0 

.55208*70 ■ 

•51752+nU 

.49716*10 

.49048*1(1 


.77319*00 

.70643*na 

.77il9*-10 

.71764*00 

.65528*00 

;6n228*OU 

•56427*00 

.54188*00 

.53452*00 


.83893*00 
.R6429+0C 
.8 309 3*C0 
.77834*00 
.7 103C*0U 
.65249*00 
.61102* )0 
.58660*00 
.57857*00 


. 90467*00 
.93214*00 
.90467*00 
•83903*00 
• 76533*00 
7D27D*0tr 
.65777*00 
.63131*00 
•62262*00 


15 

;ior»o*ci 

• 666 6 7 * 00 

‘.75791*00 

•82035*00 

.89973*00 

.97042+00 

.lD0nn*01 

.97042+CU 

•89973*00 
. 8203*^*00 
•.75291 + DO 
.70453*00 
.67603*00 
. 66667 * 00 



PEP^in CtiPVEO CtWTERt-Mrt- TEJT-€*^t fE*tPH--ttvy-- 


- WE Bj- W T» 


» > _a.XiL 

f n UL K JO 


LEVEL I IHeT/1 ***** 


— 

■ rs 

1 

■ ;tioco 

2 

.71<IJ-0X 

3 

.T4 77100 

4 

;2143*C0 

5 

•TB-STIDO 

;3l77*00 

lY 

15 

Y 

.1G00C*01 

lao.ao 

180.00 

180.00 

18 0.00 

180.00 

180.00 


"^77857♦0t3■ 

787 . 79 

18 7 . 74* 

lb 7. 14 

16 7 • 74 

7672-74 

76"7” .14 

13 

.05714*03 

154.29 

154.29 

154.29 

154.29 

154.29 

154.29 

12 

.7857 1*0C 

14 1.43 

141.43 

141.43 

14 1.43 

141.53 

141.43 

11 

.71429*03 

128.57 

128.57 

128.57 

128.57 

128.57 

128.57 

I D 

.642P6*on 

115.71- • 

115.71 

115.71 

115.71 

115.71 

115.71 

. 9 

.5714 3*00 

1Q2.R6 

102.86 

102.86 

102.86 

102.86 

102.86 

8 

•5000C*CO 

9E.OO 

90.00 

90.00 

90.00 

9C.00 

90.00 

7 

.42857*00 

77.14 

77.14 

77.14 _ 

77.14 

77. 14 

77.14 

b 

i J57 14*UU ' 

64 .77 

”■ 64 ;29 


64“.29 ■■ 

6-4279 — 

64 . 29 - 

5 

.2857 i*nn 

51.43 

51.43 

51.43 

51.43 

51.43 

51.43 

H 

.21429*00 

38.57 

38.57 

38.57 

38 .57 

38.57 

38.57 

3 

.14206*QG 

25.71 

25.71 

25.71 

25.71 

25.71 

25.71 

2 

.71429-01 

12.56 

12.88 

-72.86 ' 

12.86 

' 12.86 

12 .86 

1 

.03000 

.ECl 

.00 

.00 

.00 

.00 

.00 


7 8 lO 

.IZBIPDTI ■■lEOtnJPTra V5TlT*mJ m7V*OT 


180.00 

180.00 

180.00 

180.00 

lb/ *14 


lb f • 14 

*l6 / • X4 

IBM *29 

154.29 

154.29 

154.29 

141.43 

141.43 

141.43 

Im.43 

128.57 

128.57 

128.57 

128.57 

115.71 

- 115.77 - 

- il5.71 

115.71 

102.86 

102.66 

102.86 

102.86 

90.00 

90.00 

90.00 

90.00 

77. IM 

77.14 

77.14 

77.14 

~“64 229-- - 

tn,Z9 — 

79 

6Ti2^- 

51.43 

51.43 

51.43 

51.41 

38.57 

38.57 

38.57 

38.57 

25.71 

25.71 

25.71 

25.71 

12.86 

■ 12.86 " 

- 72.86 - 

- 12i8S 

»Q0 

.00 

.00 

.00 


17: 

7 T 

15 .lDOOC+01 
1* .9285 7f07 

13 .85 7 18*0E 

I? — ;7P57PT)C 
II .7l429»Q0 


10 .84286«00 

9 ,57m3»0'] 

8 .sooncioc 

7 .42b57*U1' 

6 .357 1'4*nc 

5 ■ .2857 UO'l 

2 .71429-01 

1 .OIOOC 


11 

714J+0D • 

12 

.7857*00 

13 

.8571*00 

i428B*'TO 

.iSoo*oi 

180.00 

180.00 

180.00 

180. UO 

180.00 

lb7.14 

167.14 

167. 14 

167.14 

167. 14 

154,29 

154 .29 

154.29 

154 .29 

154.29 

141. 4T 

■ 141-.43" 

- 141.43 — 
128.57 

- ‘141,43 — 

■~i41.43 

128.57 

128.57 

128 .57 

128.57 

115.71 

115.71 

11 5. 71 

115.71 

115.71 

102.86 

102.86 

102.86 

102.86 

102.86 

90. CO 

9 0.00 

90.00 

90.00 

9^.QD 

U:\^ 

77.14 

77.14 

77.14 

77,14 

64.29 

64.29 

64.29 

64*29 

51.43 

51.43 

51.43 

51.43 

51,43 

"38 .57" 

38.57- 

’38*57 

— 38.57 

" ■ 38.57 

25.71 

25.71 

25. 7X 

25.71 

25.71 

12.56 

12.36 

12*86 

12.86 

1?*86 

.00 

.00 

• ao 

.00 

*Cu 




-pcps-i-e -e-frst — 


Unit Ul lUrO 


BWO = Of 
IDUCT : 1 


J3-TH 5T»TION , S: .40C16805*Ql 

LOCITION OF CENTERLINE '( .QaOOQQOO , .68 36 724 G-01 , . 40QlJDDDa* 01 ) OX- .50000000*00 

KOORD- I 

1FRC3 : 1 nCNT = 0 VELOCITT BOUNDARY CONDITIONS 6 6 12 3 

INACH = 1 JFRCT : 1 DENSITY BOUNDARY CONDITIONS b 11 12 20 

vELOcn T vL ti O R ui5 P EA T E D~TmroTmmTiONia coordinates ^ 

I 


level 1 U'VEL »♦*♦* 



.... 

I ■ 

2 

T ‘ " 

' "4 ‘ ‘ 

5 


7' 


— r 

- "TO - ■ 

ir 

Z = 
y 

.OUOO 

.7143-01 

. 1429 + 00 

.2143.00 

.2857*00 

,3571*00 

.4286*00 

.5000*00 

,5714*00 

.6429*00 3 


iI?859*oJ 

• ¥?.44 t*oa 
.99520+10 

*99T62*00 

.99614+00 

!9989%SS 

• IduZB +ul 
.10023+Ul 

•10043+01 

TIW5T?Xn 

.10047+01 

»lu022+01 

•10017+01 



I?S?2itoo 

13 

.857141-00 

•9974D+1U 

.99770+10 

.99860+00 

•10036+01 

:mnsi 

;iSSI?:S} 

.99993+00 

>•98991+00 

.96012+00 

•96806*00 

.89518*00 

.90993*00 

12 

.78571*00 

.IU022*01 

.10018*01 

.17007*01 

.10013*01 

.10005*01 

,99262*00 
.10005*01 - 

•• 11 

.71474*00 

.10105*01 

'.10099*01 

.10087*01 

.10078*01 

.10080*01 

.10078*01 

.10063*01 

.98226*00 

.93378*00 

10 

9 

.64286+00 
.57343 + 015 

• 102J8 + G1 
.10258+01 

.10ZQ7+01 

•10265+01 

.10204.01 

.il283+01 

.•}8^8’^*8} 

:i8il2:g} 


.10170*01 

.10262*01 

li8^3B*81 

:IBS^?:81 

•96604+QQ 

.99703+00 

8 

•sonoc+Qo 

.10187*01 

.10197*01 

.10727*01 

.10247*01 

.10259*01 

.10265*01 

.10267*01 

.10263*01 

.10234*01 
; 1024 7 *01- 


■ T" 

.47857*00 

■.ini3D*Dl‘ 

■■ .10143F01 

■.iii82*ai 

.10266+01 

.io2i5*or " 

•10303+01 

.10239*01 

'.10257*01 ' 

.10269*01 - 

■.■1D272»01 

•10331+01 

il0133*0l 

6 

.35714+00 

.10193+11 

•10211+01 

•10326+01 

.10339+01 

.10344+151 

.10265*01 

•10034+01 

5 

4 

’Aiiim 

.10309+01 

.10410+11 

.10331^01 

.lD429*0l 

.10196*01 

.10486*01 

.10429*51 

.10509*01 

.10443*01 

.10513*01 

.l0445*0l 

.10506*01 

.10435 *01 
.10481*01 

.10390*01 

.10405*01 

:i8fSi:8i 

.98643*00 

.96861*00 

3 

.14286+nO 

. 1L476+ 11 

.■n487 + ni 

.11522+01 

;iD52B+Cl 

. 1C52«*D1 

.10512*01 

•10478 fOl 

.10378*01 

.10105*01 

.9S380*D0 

i 


.10510*01 
. 10521*01 

.10514*01 

.10573*01 

.ln526+0l 

.10527+01 

.10518*01 

.10514*01 

.10507*01 

.I050?*0l 

.10492*01 

.10486*01 

.lC454.nj 

.10446*61 


;J88H§:8i 

.94407*00 

.94084*00 


— 

- T? — - 

'll 

~ 12 

n * • 

“14 

15 


4 - 

.7143+00 

. 7857+00 

• 85 71 *U0 

.9286*00 

.1001 

IS 

•IQOUQ+01 

.76862*10 

*60324+00 

.39294+00 

•16572*CO 

.ooouo 

14 

.92857*00 

.77237+10 

•6075U+00 

,39678+00 

.16783*00 

.ouooo 


■ .78J58*ia ' 

’i62225*nO ■ 

,40829*00 

.■17416*00“- 

“. ooour 

12 

.7857 1*00 

.80381+10 

.64225+00 

••42706*00 

.18363*aO 

.00000 

11 

,71429*00 

.63566+10 

.67596+00 

.45427*00 

.19581*00 

.00000 

10 

.64 266*00 

•68135+ 

.72612+10 

,49478*00 

.2124 3*oO 

.OCOOF 

' 7 

“".5714 3*0i7 

• 93‘<53^+ 1 J 

. 79 354 *30 

.55721+10 

.2426 3.30' 

.onop 

6 

.50000*00 

.9b933+C'J 

.85015*00 

.62762*10 

.3C199.U0 

.CLO')" 

7 

• 42857.09 

.96832+ -10 

.8'‘5ll*G0 

.6551 3+00 

.35950*U0 

.OUOOC 

6 

.357 14.00 

•94132+ru 

482189+00 

.63899+00 

.38346*03 

.00000 

5 ’ 

“*28571*00 

.90823*10 

.78193*00 

.60860^00 

.38171*00 

.00000 

4 

.21429*00 

•88009+10 

.74972+00 

.58 156*00 

.37089»UP 

.00000 

3 

, 14286.00 

.65994+''3 

.72604+00 

.56303*00 

.36D48*UQ 

.00000 

2 

.71429-01 

.848'JC*nO 

.7159J+OC 

.55281*00 

.35389*00 

.ocooo 

1 

. )"noc 

. «^443!+ 1j 

.71187*10 

.5494 1+iC 

.35170*30 

.Conor 



00 


^,-PEPSIe -CtWVEO CEWTEfftINE T EST-e^se- P»tPtE- t€VT — 


.92857*00 
.85714+00 
.7857 1*00 

■ .714 Z9fi]rt 
.64286*00 
.57 14 3*00 
.50000 + 00 

~;4Z8 5'7*TiTr 
.35714*00 
.2857 1*00 
.21 429*0(1 
.14286*00 
. 714 29-0 1 
.00000 


.32672-03 
.63712-03 
.91570-03 
. 1146 5-02 
. 13239-02 
.14330-02 
.14704-02 
'■.I4340*02“' 
.13256-02 
.11506-02 
.91780-03 
.63881-03 
. 32765-03 
.00000 


*r56tI9T=lV-' 
-. 14207-03 
-.17617-03 
.6098 3-04 
.77795-03 
..21 378-02 
.37940-02 
,44664-02 
'■;"35 7n6-02 ■ 
.20355-02 
,79499-03 
,89598-04 
*,18010-03 
-.15488-03 
-.22623-19 


-.64091-03 
-.10861-02 
-.1084 2-02 
-.3 3177*03 
.13750-02 
.34795-02 
.43906-02 
.3 2961*TT2“ 
.13538-02 
-.25237-03 
-.10719-02 
-.11485-02 
-.70241-03 
-.28410-19 


-.93720-03 
-.16223-02 
-.17829-02 
*.10137-02 
.85360-03 
.31809-02 
.42062-02- 
■ V3036Z-D2- 
.86714-03 
-.90852-03 
-.17653-02 

-.24226-19 


-.11328-02 
-.19938-02 
-.22360-02 
*.15099*02 
.42217-03 
.28668-02 
.39635-02 
^.'27675=02 
.4699 1-03 
-.13959-02 
-.22380-02 
-.20852-02 
-.12142-02 
-.45884-19 


{4^ - 

~T3~;8*7mno - 

12 .7857 1*00 - 

11 .71429*00 - 

10 .64286*00 - 

■ 9 '.57 14 3*0? - 

8 .5?00C*Dr; 

7 .42857*00 - 

6 .35714*00 - 

5 — ;28571*cm * 

4 .21429*00 - 

3 .14286*00 - 

2 .71429-01 - 

J .00000 


11 

.7143*00 
.70449-20 
.16326-02 - 
.309 71- "2~ 
.41026-02 - 
.42385-02 - 
.31177-02 - 
.95324-03 ■ 

. 28864-03 - 
.87 189--13 ■ 
.30742-02 • 
.44289-02 
.43831-02 ■ 
.33679-02 ■ 
. 1 7956-02 ■ 
.57815-19 ■ 


" 12 

. 7857*00 
.19353-20 
.16515-02 
i321I7**02-‘ 
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Figure 3. Transformation as an Embedding 
into Three Dimensional Euclidian Space 
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Figure 4c. Tube— Like Surfaces of Constant Pseudo— Radius 
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Figure 6. Velocity Profiles in the Entrance Region of a Circular Pipe 
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Figure 9. S— Shaped Bend with Circular Cross Section 

(FLOW LEFT TO RIGHT) 
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Figure 10a. — Cross Flow Velocity Vector in a Curved Circular Duct 
,Sta. No. 5 
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Figure 10b. — Cross Flow Velocity Vector in a Curved Circular Duct 
, Sta. No. 6 
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Figure 10c. — Cross Flow Velocity Vector in a Curved Circular Duct 
Sta. No. 7 
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Figure 10d. Cross Flow Velocity Vector in a Curved Circular Duct 
, Sta. No. 8 
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Figure lOe. — Cross Flow Velocity Vector in a Curved Circular Duct 
, Sta. IMo. 9 


78 - 01 - 128-17 


00 

SCALE! WAX • 10*-e3 

STATION NO. le READ AT ( .0000 .1270-01 .2S00401) 


ft A 



Figure lOf. — Cross Flow Velocity Vector in a Curved Circular Duct 
, Sta. IMo. 10 
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Figure lOg. — Cross Flow Velocity Vector in a Curved Circular Duct 
, Sta. No. 11 
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Figure lOh. - Cross Flow Velocity Vector in a Curved Circular Duct 
, Sta. No. 12 
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Figure lOi. — Cross Flow Velocity Vector in a Curved Circular Duct 
, Sta. No. 13 
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Figure 1 la. — Cross Flow Velocity in a Curved Elliptic Duct, Shape = 0.667 
,Sta. No. 5 
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Figure 11b. — Cross Flow Velocity in a Elliptic Duct, Shape = 0.667 
. Sta. IMo. 6 
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Figure 11c. — Cross Flow Velocity in a Curved Elliptic Duct, Shape = 0.667 
, Sta. No. 7 
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Figure 1 1d. — Cross Flow Velocity in a Curved Elliptic Duct, Shape = 0.667 
, Sta. No. 8 
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Figure lie. — Cross Flow Velocity in a Curved Elliptic Duct, Shape = 0.667 
, Sta. No. 9 
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Figure 1 1f. — Cross Flow Velocity in a Curved Elliptic Duct, Shape = 0.667 
, Sta. No. 10 
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Figure 11g. - Cross Flow Velocity in a Curved Elliptic Duct, Shape = 0.667 
, Sta. No. 11 
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Figure 1 1h. — Cross Flow Velocity in a Curved Elliptic Duct, Shape = 0.667 
, Sta. No. 12 
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Figure 1 1i. — Cross Flow Velocity in a Curved Elliptic Duct, Shape = 0.667 
, Sta. No. 13 
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Figure 12a. — Cross Flow Velocity in a Curved Elliptic Duct, Shape =1.5 
, Sta. No. 5 
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Figure 12b. — Cross Flow Velocity in a Curved Elliptic Duct, Shape =1.5 
, Sta. No. 6 
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Figure 12c. Cross Flow Velocity in a Curved Elliptic Duct, Shape = 1.5 
, Sta. No. 7 
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Figure 12d. — Cross Flow Velocity in a Curved Elliptic Duct, Shape =1.5 
, Sta. No. 8 
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Figure 12e. — Cross Flow Velocity in a Curved Elliptic Duct, Shape =1.5 
, Sta. No. 9 
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Figure 12f. — Cross Flow Velocity in a Curved Elliptic Duct, Shape =1.5 
, Sta. No. 10 
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Figure 12g. — Cross Flow Velocity in a Curved Elliptic Duct, Shape = 1.5 
, Sta. No. 11 
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Figure 12h. — Cross Flow Velocity in a Curved Elliptic Duct, Shape =1.5 
, Sta. No. 12 
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Figure 12i. — Cross Flow Velocity in a Curved Elliptic Duct, Shape = 1.5 
, Sta. No, 13 
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Figure 13a. — Cross Flow Velocity Vector in a Curved Superelliptic Duct, Exponent = 10 
, Sta. No. 5 
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Figure 13b. — Cross Flow Velocity Vector in a Curved Superelliptic Duct, Exponent = 10 
, Sta. No. 6 
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Figure 13c. — Cross Flow Velocity Vector in a Curved Superelliptic Duct, Exponent = 10 
, Sta. No. 7 
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Figure 13d. - Cross Flow Velocity Vector in a Curved Superelliptic Duct, Exponent = 10 
, Sta. No. 8 
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Figure 13e. — Cross Flow Velocity Vector in a Curved Superelliptic Duct, Exponent = 10 
, Sta. No. 9 
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Figure 13f. — Cross Flow Velocity Vector in a Curved Superelliptic Duct, Exponent = 10 
, Sta. No. 10 
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Figure 13g. - Cross Flow Velocity Vector in a Curved Superelliptic Duct, Exponent = 10 
, Sta. No. 11 
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Figure 13h. — Cross Flow Velocity Vector in a Curved Supereltiptic Duct, Exponent = 10 
, Sta. No. 12 
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Figure 13i. — Cross Flow Velocity Vector in a Curved Superelliptic Duct, Exponent = 10 
, Sta. No. 13 



Figure 14. Transition Duct Superellipse Exponent - 10 to Circular 
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Figure 15. Transition from Superelliptic Pipe (C = 10) to Circular Pipe 
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